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ABSTRACT 


This document describes the parameterization of longwave radiation in the 
UCLA/GLAS general circulation model. Transmit tances have been computed from 
the work of Arking and Chou for water vapor and carbon dioxide and ozone 
absorptances are computed using a formula due to Rodgers. Cloudiness has been 
introduced into the code in a manner in which fractional cover and random or 
maximal overlap can be accomodated. The entire code has been written in a form 
that is amenable to vectorization on CYBER and CRAY computers. Sample clear 
sky computations for five standard profiles using the 15- and 9-level versions 
of the model have been included. 



1. INTRODUCTION 


The parameterization of longwave radiation for general circulation models 
is currently under extensive study. Several experiments have shown that model 
integrations over long periods are very sensitive to the particular methods 
used in calculating the radiative terms. However, because of the spectral 
nature of radiation, an accurate description of the radiation field requires a 
summation over several spectral bands. Moreover, since the Interaction paths are 
very long, transmittances between each atmospheric level and all other levels 
must be computed and carried through the spectral integration. Because of 
this, the computation of longwave radiation is the most time consuming step in 
the Integration of a general circulation model. One common practice Is to use 
a fairly detailed and accurate representation of the radiation, but to update the 
model infrequently, say every five hours, even If the physics is updated every 
half hour or hourly. The problem with this approach is that the life-cycle of 
clouds is often less than the time-interval of the update, and so there can be 
no interaction between radiation and cloudiness on shorter time scales. The 
need, therefore, is for a fast yet accurate treatment of longwave radiation. 



2. RADIATIVE TRANSFER 

For a clear atmosphere, the upward and downward fluxes at pressure level p, 
integrated over the spectral range Av, can be written as 

P dx (p,pO 

F clr + ( P> = J dv [®v < T S ) T v ( P*P S ) + / B v < T (pO) — dp'l 

Av p s dp (la) 

and 


, r P d t (p ") 

F clr + <P> = / dv I / B v (T( P'> —T? d P'l 

Av P+. p 


(lb) 


where B V (T) is the blackbody flux at temperature T and wavenumber v, x v (p,p") 
is the diffuse t ansmittance between the levels p and p", and the subscripts s 
and t denote the surface and the top of the atmosphere, respectively. In order 
to use Equations (la, b) as written, the spectral interval must be very narrow, 
since t v varies very rapidly with wavenumber. In practice, models avoid this 
problem by resorting to certain assumptions regarding the distribution of spectral 
lines in the various absorption bands and constructing approximate band models. 

We have adopted the method described by Chou and Arking (1980) and Chou (1984) 
for the water vapor bands, Chou and Peng (1983) for carbon dioxide and Rodgers 
(1968) for ozone. Table 1 lists the spectral regions covered for each of the 
bands in the model. Apart from line absorption in the various bands, there is 
also water vapor continuum absorption over a wide spectral region that overlaps 
with water vapor line absorption as well as the carbon dioxide and ozone absorption 
bands. This continuum absorptance, computed from a formula given by Roberts 
et al . , (1976) Is multiplied by the line transmittances where appropriate. 

Details are given in Chou and Kouvaris (1981). 

An alternate form of the flux equations is obtained through integration by 
parts of (la,b) reducing them to 
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T(P S > 

F c l r + (P) = B(T(p)) + G (p,p s ,T s ) - G (p,p s ,T(p s )) + J 3G(p,p',T(p ') ) / 3T dT(p') 

T(p) 


(2a) 

T(p t ) 

F cir + (p) = B(T(p) ) - G (p,P t ,T( Pt )) + / 3G(p ,p %T(p ') )/ 3T dT(p') , (2b) 

T(p) 


where 

G (p,p',T) - / t v (p,pO B v (T)d v , (2c) 

Av 


3G(p,p',T)/3T = / T v (p,p') 3B V (T)/3T dv 

Av 


(2d) 


In the 9,6 un and 15 vim bands, the spectral width of the band is sufficiently 
narrow such that a mean value of the Planck function may be used. That is, 

G(p,p',T) = / t v (p,p')dv/Av •/ B v (T)d v (3) 

Av Av 

For carbon dioxide, the spectral integration is performed separately over the 
band center from 620-720 cm'^ and the band wings from 540-620 cm~l and 720-800 cm~l. 
The transmittance in the ozone band is obtained by scaling the band absorptanca 
given by Rodgers (1968), assuming that all absorption by ozone is confined to 
the 980-1100 cm~^ region. For water vapor, the variation of B v and t v in (2c) 
is accounted for over the various spectral ranges. Details of the technique 
used to compute transmittances may be found in the references quoted above. 

For cloudy skies, If an assumption is made that the clouds are black at 
thermal wavelengths, the above equations can be written very easily for a 
completely overcast sky by replacing the limits on the exchange Integral, since 
transmission through a black cloud is zero. However, we have retained the 
possibility of allowing fractional cloud cover In the various model layers or 


- 3 " 



partial transmittance through high tenuous clouds. Therefore, a slightly 
different computational approach has been used utilizing the probability of a 
clear line of sight between each layer and all other layers, the ground and the 
top of the atmosphere. This will be explained in the description of the code. 

The final result required by the dynamic model is the radiative cooling 
rate, or the divergence of the net flux. In terms of the quantities computed by 
Eqs. (2a, b) the cooling rate is 

- dT = _£ _d ( F+ _ Ff ) (4) 

dt Cp dp 

Of course, all other fluxes are available and the net flux at the surface is 
needed to perform the surface energy balance. For diagnostic purposes, the 
upward flux at the top of the atmosphere is also retained as output of the 
radiation model. 
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3. RADIATION CODE 


In principle, if the atmospheric parameters such as temperature and mixing 
ratio of water vapor, carbon dioxide and ozone are given, the various transmit- 
tances required for the flux computation can be obtained using the theory 
covered in the various references cited earlier. In general, the transmittance 
depends on the atmospheric temperature, presure and amount of absorbing gas. 

The atmosphere, however, is inhomogeneous in that the absorber amount, pressure 
and temperature are functions of altitude, and the absorbing path has to be 
suitably scaled in order to compute the transmittan. e . Again, the references 
cited explain the theory and techniques utilized. Practical difficulties arise 
when the exchange integrals in (2) are to be set in a numerical procedure to 
compute fluxes. The dynamic atmospheric model provides only the surface temper- 
ature and layer mean temperatures and mixing ratios for the model layers. However, 
the flux integrals are over a continuous atmospheric path, and the modeler is 
faced with the familiar problem of discretization. Since transmittance is a 
non-linear function of atmospheric path, widely varying results can be obtained 
depending on the assumptions made regarding the structure of temperature and 
mixing ratio between model levels. This problem of interpolation h^s been dealt 
with In various ways. Usually, the dynamic model layers are further subdivided 
and assumptions made about the profiles within the layer. Even after this, 
different quadrature schemes can yield different results, especially over paths 
where the transmittance varies rapidly. 

In the radiation code that we have designed, we have side-stepped the inter- 
polation problem by computing absorber amounts assuming that the mixing ratio 
of absorbing gas In each atmospheric layer is constant. The flux computations 
require that the temperature be known at the layer interface (each level of the 
model). For this, we interpolate the tempers ture linearly In the logarithm of 



pressure for all levels except the top and bottom of the planetary boundary 
layer which is the lowest layer of the model. The temperature at the top and 
bottom of this layer is computed by following the appropriate adiabat. The 
dynamic model furnishes information up to the highest layer in the model, but 
there is some residual opacity in the upper atmosphere that, if ignored, could 
influence fluxes and cooling rates in the upper layers of the model. Usually, 
some radiation layers are inserted above the model layers and climatological 
information regarding temperature and mixing ratios supplied to the radiatxc 1 
code. In the 15-level UCLA model, the uppermost layer is centered at 1.4 mb 
with its bottom at 1.9 mb, and its upper boundary at 1 mb. We insert one 
radiation layer between 0 and 1 mb, and specify the temperature and mixing 
ratios in this layer, for example that they are the same as those in the top 
layer of the model. 

The description of the radiation code is given below in separate sections 
for each subroutine and a listing Is attached in the appendix. The code as pre- 
sented is designed for clarity rather than optimization. It is also written 
so that it can be vectorized quite easily for use on CYBER and CRAY machines. 

Only key statements are explained with references to the papers from which the 
equations were obtained. 

SUBROUTINE RADIR 

This is the driving program of the longwave radiation code. The Input to 
this code is the atmospheric profile of temperature and gaseous constituents. 

The input parameters passed to the subroutines are copied into internal variables. 
NDATA is the number of layers in the model which is 15 for the 15-level UCLA 
GCM. The other variables are as follows. 

PU = pressure at each level of the model 

PA = pressure at the middle of each layer 
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TA = mean temperature of the layer defined at pressure PA 
TS = surface temperature 

WA = water vapor mixing ratio of the layer defined at pressure PA 

OA = ozone mixing ratio of the layer defined at pressure PA 

CSS = random overlap cloud fraction 
CCU * maximum overlap cloud fraction 

The pressure at the top and bottom of the extra layer inserted on top of 
the dynamic model is designated PU(1) and PU(2) and defined. The cloud fraction 
in this layer is 0.0 and all other cloud fractions are constrained to lie between 

0.0 and 1.0. After calls to the other subroutines, the fluxes at each level 

are returned to RADIR for the computation of cooling rates using (4). The term 
FAC is a conversion factor and since fluxes are in e.g.s. units, a conversion 
is also applied to obtain the upward flux at the top of the atmosphere (ULWTOP), 
the downward flux at the surface (DLWBOT) and the net upward flux at the surface 
(RS) in Wm”2. The cooling rates in the 15 model layers are returned to the main 
program and are in °C day"*. 

SUBROUTINE INDATA 

When the program is run for the very first time, data is read in through 
this subroutine call. The data concerns the G-f unction and temperature 
corrections in the water vapor bands and is stored in tables. There is a 
subroutine call to a table look up program that computes the relevant numbers. 
This will discussed later. 

SUBROUTINE SETUP 

This routine does some housekeeping on the input variables and also performs 
the tempeiature interpolation and computes scaled absorber amounts. At the end 
of the interpolation, the variables are on a vertical grid as shown in Figure 1. 
The amount of water vapor and ozone in each layer is computed from 
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( 5 ) 


2 

U(p 1 ,p 2 ) - i / q dp - J Ap 1>2 

g Pi 8 

where q Is the mixing ratio. 

A factor of 1.02 is used to convert the pressure to c.g.s,. units and divide in 
the acceleration constant. Some points to note are that there is a lower limit 
on the allowed mixing ratio. This is to prevent numerical problems at later 
stages . 

Although (5) gives the amount of water vapor in a column of air from 
p to p^ as mentioned earlier, it is necessary to scale this amount to compute 
the transmission. In the routine we have defined three scaled water vapor 
absorption path:, one each for the band centers and band wings and one for 
e-type absorption. The path for band abso. ^tion is of the form 

(Pi,P 2 ) - f ^ (P"/P r )“ R [T(p"),T r ] q(p") dp*7g ( 6 ) 

P' 

and for e-type absorption is 

P 

w (p ,p ) - / e(p") . exp [1800 (1/T(p") - 1/T 0 )1 q(p”) dp’Vg. (7) 

3 12 p - 

The temperature dependent function R in ( 6 ) is written as SZ1 for the band 
center and SZ2 for the band wings and is given by 

R(T,T r ) - exp [r(T - T r )] , . 8 ) 

where r ■ 0.005, T r ■ 225°K for the band centers and r * 0.016 and T r * 256 # K 

for the wings. The reference pressure, p r , is 275 mb for the center and 550 mb 

for the wings and m * 1 for both. In (7), e(p ) is the water vapor partial 

pressure in atmospheres and T 0 * 29C*°K. The equation is coded as XX * ........ 

with UBAR denoting the scaled path length, W 3 , from each pressure level to the top 
of the atmosphere. The factor 630 in the code converts mixing ratio to vapor pressure 

in atmospheres and is (0.622*1013) and ex , j lined in Chou and Kouvaris (1981). 
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VBAR and WBAR are respectively, the scaled path le-.gths in the band center and 
wings. UBARM, VBARM and WBARM are the means of path lengths from two adjacent 
layers to the top, a quantity that is required later in the computation of 
transmission. 

The scaled amount of carbon dioxide denoted by PSCALV for the band center 
and PSCALW for the wings is obtained from (6) with c(p") replacing q(p") where 
the concentration of CO 2 is in ppm. A conversion factor of 0.26 for 330 ppm is 
used here. Any other concentration can be used by changing the number accordingly 
since scaling in c is linear. The reference pressure, p r , is 30 mb in the center 
and 300 mb in the wings while m = 0.85 in the center and 0.5 in the wings. For 
pressures p<p c , the pressure scaling is (p c /p r ) where p c * 1 mb. The temperature 
scaling R(T,T r ) is written as in (8) with r = 0.0089 in the center and 0.025 
in the wings; T r = 240°K for both regions. As in water vapor, the mean scaled 
path lengths PSCAMV and PSCAMW are also computed. 

The ozone amount is given by (5) without any temperature scaling replacing 
q by the ozone mixing ratio. Additionally a quantity called the effective 
broadening pressure for a path is required and this is given by 



q p dp// 


q dp 


(9) 


The ozone amount is coded as V while the numerator of (9) is coded as U in the 
subroutine. Additionally mean quantities UM and VM are also computed. 


SUBROUTINE IRFLUX 

This subroutine computes the fluxes in each of the bands in turn and adds 
them up. The various terras in (2a,b) are computed by calls to other subroutines 
that calculate the appropriate transmission functions. 

The G-f unctions for water vapor defined by (2c,d) are computed here* As 
explained in Chou (1984), 
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( 10 ) 


G(w ,w ;T) - T (w ,w ;T) B (T) 

13 13 

where wj, W3 are defined by (6) and (7) and B(T) is the Planck function at T 

integrated over the wavelength interval in question. The transmission, t(w ,w ;T) 

1 3 

is itself parameterized as 

t ( w i» w 3» t ) - t (w 1 ,w 3 ; 250)[1 + a (wj,w 3 )(T - 250) + 8 ( w i» w 3 )(T - 250) 2 ] 

(ID 

The quantities, x (w ,w ; 250), c and 8 are obtained from a table lookuo through 

1 3 

a call to subroutine TABLE or a parameterization in subroutine GFUNC. For the 
band cencer, only wj is needed whereas for the band wings both wj and W3 are 
needed. Note that wj is scaled differently for the center and wings, such that 
the passing argument to TABLE or GFUNC contains chree path lengths UBAR,VBAR 
and WEAR. 

The downward flux is computed from (2b) for which the terms G(p,p t ,T(p t )) 

Pt 

and l AG(p,p ,T(p ))) are required. In the code, G(p,p t ,T(p t ) is called SHT 
P 

and is computed for each level in the model starting with the second level and 
proceeding down to the surface. 3y definition, the downward flux at the topmost 
radiation level, 0 mb, is zero and so no computations are performed for that 
level. For each level, the term AG is computed with respect to all other layers 
following an approximation given in Chou and Kouvaris (1981). Basically, the 
approximation states that the contribution to the exchange integral at level I 
from a layer J is the product of the change in Planck function across the layer 
J and the transmission from the middle of layer J to level I hence the need for 
computations of mean absorber amounts in SETUP. The terms AG are called SG in 
the code for both upward and downward fluxes. Note that SG is two-dimensional 
and only one half of the array is used for the downward flux, whereas the other 

half is used for the upward flux. For the upward flux, [G(p,p s ,T s ) - G(p,p 8 ,T(p s ))] 

is called SHU and is computed for each level except the surface. 
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The Planck function integrated over the spectral range in each band is 
computed by linear interpolation from tables given by DATA statements* The 
nomenclature is 

SV * Planck function in water vapor band center 

Stf * Planck function in water vapor band wings 

BAI * Planck function in 15 pm region 

BWI = Planck function in 9.6 m region 

SH = Planck function integrated over entire spectrum 

The terms in the flux equation in the 9.6 and 15 m region are competed 
after calls to subroutine TRC03 which returns the transmission in the 15 un 
region as TAUF and in the 9.6 \*n region as TXUF. These, when multiplied by the 
appropriate difference in Planck function across the layer are added to the SG 
calculated for water vapor. 

After G and AG are computed, the upward and downward fluxes, FU and FD are 
calculated. For the clear sky, all cloud parameters CS, CC and CH are unity and 
the fluxes are simply summations of AG over the appropriate layers. The variable 
SG is summed in SUM and the aownward flux is written as 

FD = SUM - SHT + SH 
while the upward flux is written as 

FU = SUM + SHU + SH 

following (2a,b). The upward flux at the surface is the black body flux since 
the surface is assumed to have unit emissivity in this code. Therefore, for 
the surface, FU - SH. 

SUBROUTINE TRC03 

This subroutine computes the transmittance in the CO 2 and O 3 regions given 
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the scaled atmospheric paths for the various absorbers. 

The transmittance in the 15 m region is expressed as 

T (Hl.W 2 .W 3 ) = T (wj) • X(W2> * T (W 3 ) (12) 

where subscripts 1 , 2 and 3 refer to water vapor line, CO 2 line, and water 
vapor e-type absorption, respectively and the multiplicative property of 
transmittances has been used. Expressions for t (w^) and x (^ 3 ) in this region 
are given in Chou (1984) and are 

x (w^ = exp [-6.7^/il + (13) 

and 

t (w 3 ) = exp (-27 w 3 0 * 83 ) (14) 

The CO 2 band transmittance is expressed as 

x (w) = exp [-aw/(l + bw n )] (15) 

for both the band center and wings. The values of a, b aud n are given in Chou 
and Peng (1983). Note that the wide band model has been used here. The trans- 
mittance x (W 2 ) is given by suitably weighting (15) in the band center and 
wings, i.e. 

T (w 2 ) = 1- [t c Av c + -^ Av w ] (16) 

where c and w refer to the center and wings respectively. For the model used here, 
Av c = 100 cm 1 and Av w = 160 cm ^ and the ratio of weights is therefore 5 : 8 . 

The transmittance is denoted by TAUF and passed back to IRFLUX. 

The transmittance in the 9.6un region is computed by a modification of the 
band absorptance formulation given in Rodgers (1968). If we assume that all 
the ozone absorption is restricted to the 980-1100 cm“l interval, then the 
average transmittance in that interval based on Rodgers 1 formula is 
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1 - PI - 2 . 1 [1 - exp (-4.398 p 
120.0 


. / 345.28U 0 , . 

{/l + 4 ( 23 ) - 1} ] 


0.876? 


( 17 ) 


-2 


where Uq^ is the column amount of ozone expressed in gem and p is the effective 
broadening pressure in atmospheres. If p is the mean pressure of the atmospheric 
layer and Ap is the pressure thickness of the layer in mb and X is the ozone mixing 
ratio in g/g in the layer, then 


U 0 ^ - 1-02 X Ap 


(18) 


for each layer and the path through several contiguous layers is 

U 0 = l 1.02 X* Api . (19) 

3 i 

The effective broadening pressure in atmospheres, for a path through several 
contiguous layers is 

p * [^ X i P i Ap i]/iQi3.25 (20) 

Wi Api 

The ozone amounts and the numerator of (20) were calculated in SETUP and are 

passed from IRFLUX. A note of caution here is that if 0ZAI in SETUP is very 

small, there is a danger of computational overflow in using (20). This can be 

rectified by setting a small non-zerc value for OZAI if OZAI is zero or very 

nearly so. If the model is to be run without ozone, XZ * 1.0 in the subroutine. 

The transmittance in the ozone band has to be multiplied by the transmittance 

in the water vapor continuum that overlaps this region. The weak water vapor line 

absorption it neglected so only w computed by (7) is needed. The transmittance 

3 

is given in Chou (1984) and is 

t(w ) - exp (-9.79v ) (21) 

3 3 

and is denoted by TXUF. 
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SUBROUTINE TABLE 


Calls to this subroutine are aade to compute t(w ,v ; 250 ), a(w ,w ) 

13 13 

and 8(w ,w ) as defined by (11). Tables are read in at the start of the 
1 3 

program and during an integration the values are interpolated and returned to 

the calling subroutine. A description of the table and the tables themselves 

are given in Chou (1984). The table is used ouly for the water vapor bands. 

For the band center, the quantities ~re a function only of VBAR which is w for 

the center, called B in the subroutine, t, a and 8 are given at regular Intervals 

in log w ; there are 22 values from log w » -6.0 to log w * 0.3. Extra- 
10 1 10 1 10 1 

polation is permitted for larger values of log w but a check is made to ensure 

10 1 

that T 0.0. t, a, 8 are returned as WV, YV, ZV, respectively. 

For the band wings, both w and w determine x, a and 8 and so the table is 

1 3 

two-dimensional with 22 values of log w from log w = -5.4 to log w =0.9 

10 1 10 1 10 1 

and 19 values of w from 0.0 to 0.054. The scale for w is linear unlike that 
3 3 

for w . In the wings, w Is WBAR, called A in the subroutine, while w is VBAR, 

1 l 3 

call C in the subroutine. Again a check is made to ensure that t >_ O.*' The 

two-dimensional table lookup is linear in log w and w . x, a, 8 are returned 

10 1 3 

as WW, YW, ZW, respectively. 

Since the table lookup procedure cannot be vectorized very easily, we have 
tried to replace the table by polynomial fits valid in the range of the table. 

This subroutine called GFUNC is also provided in the code and has the same input and 


output arguments as TABLE. 


SUBROUTINE CLOUDS 

In the longwave, clouds are considered to be non-reflecting and are assumed 
to completely fill a model layer in the vertical. In the horizontal, clouds 
are permitted to occupy a fraction of the grid box or to be partially transmitting. 
In either case the cloud fraction is expressed as a number between 0.0 and 1.0, 
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with 0.0 corresponding to clear skies and 1.0 to a completely filled opaque 
cloud. A completely filled layer of thin clouds of given emissivity is designated 
by the same parameter between 0.0 and 1.0 since for non-reflecting clouds, 
partial filling and partial transmission are interchangeable. The cloudy sky 
fluxes a computed from the same set of equations as the clear sky with certain 
multiplicative cloudiness factors depending on probabilities of clear line of 
sight between levels. These probabilities are computed in SUBROUTINE CLOUDS, 
the input. :o which specifies -loud fraction in each level and also whether the 
clouds are randomly overlapped or maximally overlapped. In the program the 
two /inds jf cloudiness are specified independently, with CSS referring to 
randon overlap fraction and CCU to the maximum overlap fraction. 

The flux equations for partially cloudy skies are given below for the 
simple case illustrated in Figure 2. Consider a 5 layer model with a cloud in 
layer 3, If that is the only cloudy layer and fractional cover is N, then (2b) 
for the downward flux at a level below the cloud can be written as 


T(p t ) 

F cld +<P> = BfT(p) ) - v* ' N) G(p,p t ,T(p t )) + (1 -N) / 8G(p,p',T(p'))/ 3T dT(p') 

T(p) 

r T(p c. } 

+ N / b 3G(p,p',T(p'))/3T dT(p') (22) 

T(p) 

where p c is the level of the cloud base, in this case, level 4. But (22) can 
b 

be written as 

T( j 

F cld +<p) • B < T <P - (1 - N) G(p,p t ,T(p t )) + (1 - N) / aG(p,p',T(p'))/3T dT (p') 

T(p cJ 

- ( p c ) b 

■* J b 3G(p',p',T(p'))/3T dT(p') (23) 

T(p) 

by sp itting the limits of the first integral. Noting that G(p,p t ,T(p t )) 
involves the transmittance between level p and the top of the atmosphere, p^, 
it is evident that (1 - N) is the probability of a clear line of sight between 
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p and p t since p is below cloud base. Similarly for all levels between p c 

b 

and p t , the probability between p and p' is also (1 - N). For levels between 

p and p c the probability is 1.0. Therefore (23) can be written as 
b 


F c id +<P) 


T(p t ) 

B(T(p)) - C(p,p t ) G(p,p t ,T(p t )) + / C(p.p') 

T(p) 


3G(p,pM(p'))/3T dT(p ') 
(24) 


where C(p,p") is the probability of a clear line of sight from p to p'. Comparing 
(24) with (2b), we see that the form is the same except for the multiplication 
of the factors C(p,p t ) and C(p,p'). These are coded as CH and CS respectively. 

It can be shown that the upward flux above the cloud layer is given by 


F cld *(P> = B(T(p)) + C(p,p s ) [G(p,p s ,T S ) - G(p,p s ,T( p s ) ) ] 

T(P S ) 

+ / C(p.p') 3G(p,pM(pO)/3T dT(p') (25) 

T(p) 

The probability of a clear line of sight from level p to the surface is coded 
as CC. These cloud factors CH, CC and CS are passed back to IRFLUX for flux 
computations. 

If there are clouds in more than one layer, (24) and (25) can still be 
used but the cloudiness factor will depend on the fraction in each layer and 
cloud overlap. The problem oi fractional cloud overlap is still under study 
and there is no consensus on the preferred form of overlap in general circulation 
models. This question is only of academic importance in models that assume 
either clear skies or complete overcast. In the current radiation code we 
allow for two kinds of overlap, random and maximum. These are illustrated in 
Figs. 3 and 4. 

Only some selected probabilities of clear lines of sight will be illustrated. 
In Figs. 3 and 4 for example, Nj_, N 2 and N 5 * 0.0, and C( 6 ,l) can be seen to be 
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C( 6 , 1) - MIN [(I - N 3 >, (1 - N 4 )] 


(26) 


for maximum overlap, where MIN is the minimum function. For random overlap, 

C( 6 , 1) - (1 - N 3 ) * (1 - N 4 ) (27) 

In obtaining C for downward flux, the probabilities are computed starting from 
the top going down using equations like (26) and (27). In general, cloud 
factors are of the form where i is the level at which the flux is to be 

computed and j is the level to which the probability of a clear line of sight 
is needed. For the downward flux, j = 1 corresponds to the factor CH in the 
code while for the upward flux if j Is the lowest level of the model (surface), 
then the factor is CC in the code. The downward cloud factors are computed by 


C i>:j = MIN [(1 - Nj), (1 - Nj +1 ) , (1 - Nj +2 ) (1 - Ni-i)), (28) 

for maximum overlap and 

C i}1 = (1 - Nj) * (1 - N j + 1 ) * (1 - Nj +2 ) * * (1 - N^) (29) 

for random overlap. The upward cloud factors are computed by 

Cij - MIN [(1 - N^), (1 - Nj. 2 ), (1 - N j — 3 ) , (1 - N t )j, (30) 

for maximum overlap and 

C i,j * (1 ' N j-l) * (1 - Nj_ 2 ) * (1 - Nj-3 ) * * (1 - N t ) (31) 


for random overlap. The cloud factors for maximum and random overlap are computed 
separately and then multiplied to give the final probability of clear line of 
sight. It Is possible to have a maximum overlap cloud fraction and random overlap 
cloud fraction assigned to the same layer but it Is unlikely that a model will 
be able to make such specifications. We anticipate that the most common usage 

will be for layers of maximally overlapped cloud fractions in consecutive layers 
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overlapped with a broken cloud layer or cirrus layer separated by clear layers* 
For example in Fig. 3, if was the fraction of randomly overlapped clouds in 
layer 1 with N 3 and N 4 the maximum overlap cloud fractions in layers 3 and 4 
respectively and N 2 = 0 signifying a clear layer 2 , then (26) would read 

C( 6 ,l) = MIN [(1 - N 3 ), (1 - N 4 )] * (1 - N X ) (32) 

However, if the downward flux at level 3 was needed, the appropriate factor 
would be C(3,l) and this would be given by 

C(3,l) = (1 - N X ) (33) 

The cloud fractions in layers 3 and 4 would not enter into the computations 
since level 3 is above the two layers. This simple example illustrates the 
usage of SUBROUTINE CLOUDS. 
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4. SAMPLE COMPUTATIONS 

The radiation code has been used to compute fluxes and cooling rates using 
five standard profiles shown in Tables 2-6. The input profiles are currently 
being used in an intercomparison project in which GLAS is a participant. The 
input parameters have been interpolated to the levels corresponding to the 15- 
and 9-layer versions of the UCLA/GLAS GCM. 

The momenclature used in the tables is as follows . 

PA = pressure at the middle of each layer in mb 

PU = pressure at the bottom of each layer in mb 

TA = temperature of the layer defined at pressure PA in Kelvin 
QL = water vapor mixing ratio in layer in g/g 

03 = ozone mixing ratio in layer in g/g 

UP15 = upward flux at each level in 15-layer model in Wm“2 
UP9 = as above in 9-layer model 
DN15 = downward flux at each level in 15-layer model in Wm~2 
DN9 = as above in 9-layer model 
NT15 = net upward flux at each level in 15-layer model in Wm~2 
NT9 = as above in 9-layer model 

CR15 = cooling rate ^n each layer in 15-layer model in Celsius/day 
CR9 = as above in 9-layer model 

The same profiles have been used for the two models. The 9-layer model has 
the same vertical structure as the nine lowest layers of the 15-layer model. The 
top of the 15-layer model is at 1.0 mb whereas the top of the 9-layer model is 
at 51.8 mb. Fluxes and cooling rates have been presented such that the values 
in layers or levels at the same pressure are in the same row in the table. Of 
course, results are not presented for the six upper layers missing in the 9-layer 
model. The cooling rates computed for the 9-layer model do not differ substantially 
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from the 15-layer results except for the top one or two layers of the model where 
the cooling rate itself is quite small. These rasults also compare favorably 
with results obained for 31 layers using the same profiles# The version of 
the code used to obtain the results in the five tables used SUBROUTINE TABLE 
for the water vapor transmittance and computed the temperature at the top and 
bottom of the boundary layer by following a dry adiabat# However, the code given 
in the appendix uses interpolated temperatures for these two levels. 
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TABLE 1 


Spectral regions in the longwave radiation scheme. 



H 2 O band 
centers 

H 2 O band 
wings 

CO 2 band 
center 

CO 2 band 
wings 

O 3 band 

Spectral 

range 

(cm -1 ) 

0-340 

1380-1900 

340-540 

800-980 

1100-1380 

1900-3000 

620-720 

540-620 

720-800 

980-1100 
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TABLE 2 


TROPICAL PROFILE 


L 

PA 

PU 

TA 

QL 


03 

1 

l.il 

1.9 

267.8 

0.325E 

‘05 

0.557E-05 

2 

2.7 

3-7 

258.2 

0.325E 

-05 

0.890E-05 

3 

5.2 

7-2 

2A7.6 

0.325E 

-05 

0.105E-0L 

A 

10.0 

13-9 

237.8 

0.325E 

‘05 

0.123E-0L 

5 

19.3 

26.8 

228.6 

0.325E 

-05 

0.101E-0L 

6 

37.3 

51.8 

216.3 

0.325E 

-05 

0.505E-05 

7 

72.0 

100.0 

201.4 

0.325E 

-05 

0.920E-06 

8 

125.7 

152.6 

202.0 

0.331E 

-05 

0.216E-06 

9 

187.2 

223. L 

218.3 

0.742E 

-05 

0. I49E-06 

10 

268.2 

311*. 9 

233.9 

O.786E 

-04 

0.999E-07 

11 

372.1 

L3i.it 

249-3 

0.LL3E 

-03 

0.751E-07 

12 

502.2 

575.5 

265.0 

0.143E 

-02 

O.656E-07 

13 

662.1 

751.6 

279.6 

0.398E 

-02 

0.588E-07 

1J» 

855.7 

963.0 

291.2 

0.109E 

-01 

0.561E-07 

15 

987-9 

1013.0 

298.7 

0.158E 

-01 

0.A90E-07 


SURFACE TEMPERATURE - 300.0 


FLUXES COOLING RATES 

WATTS PER SQUARE METER DEGREES / DAY 


L 

UP15 

UP9 

DN15 

DN9 

NT15 

NT9 

CR15 

CR9 

1 

297.1 


1.1 


296.0 


5-79 


2 

297.0 


1 .6 


295.4 


5.28 


3 

296.7 


2. L 


29L.2 


3.71 


k 

296.2 


3.5 


292.7 


2.37 


5 

295.9 


5.1 


290.8 


1.16 


6 

296.0 


7-0 


289.1 


0.01 


7 

297.6 

296.7 

8.5 

5.2 

289.0 

291.5 

-0.03 

O.OL 

8 

299.6 

299.6 

10.3 

8.3 

289.2 

291.3 

0.09 

0.17 

9 

302.5 

302.5 

13.8 

12.3 

288.7 

290.2 

0.L3 

0.L7 

10 

308.1 

308.1 

23.0 

21.8 

285.1 

286.3 

1.79 

1.81 

ll 

318.3 

318.3 

52.7 

51.6 

265.6 

266.7 

2.39 

2.39 

12 

336.5 

336.5 

103.8 

102.8 

232.7 

233.'/ 

2.17 

2.17 

13 

364.5 

36A.5 

168.9 

167.9 

195.6 

196.6 

2.18 

2.19 

lit 

A00.2 

LOO. 2 

250.2 

249.2 

150.0 

151.0 

2.L7 

2.L8 

15 

447.2 

LL7.2 

359.0 

358.3 

88.2 

88.9 

3.53 

3-55 

16 

1*59-3 

L59.3 

392.0 

391 .4 

67.3 

67.9 





TABLE 3 


MIDLATITUDE SUMMER PROFILE 


L 

PA 

PU 

TA 

QL 

03 

1 

1 .it 

1.9 

271*. 7 

0.A00E-05 

0.606E-05 

2 

2.7 

3.7 

26A.8 

O.AOOE-05 

0.103E-0A 

3 

5.2 

7.2 

252.5 

0.A00E-05 

O.MAE-OA 

A 

10.0 

13.9 

2A1.5 

0.A00E-05 

0.1A7E-0A 

5 

19.3 

26.8 

232.2 

0.A00E-05 

0.10AE-0A 

6 

^ 7.3 

51.8 

222.1 

0.A00E-05 

0.583E-05 

7 

72.0 

100.0 

216.9 

0.A00E-05 

0.270E-05 

6 

125.7 

152.6 

216.1 

0.400E-05 

0.955E-06 

9 

187.2 

223-A 

217-7 

0.819E-05 

0.A71E-06 

10 

268.2 

31M 

233.1 

0.113E-03 

0.2A0E-06 

11 

372.1 

A31.A 

2A8.0 

0.A02E-03 

0. 15 ii*o6 

12 

502.2 

575.5 

262. A 

0.103E-02 

0.103E-06 

13 

662.1 

751.6 

275.6 

0.292E-02 

0.756E-07 

lit 

855.7 

963.0 

287.8 

0.731E-02 

0.582E-07 

15 

987.9 

1013.0 

293.1 

0.110E-01 

0.515E-07 


SURFACE TEMPERATURE - 29A.0 


FLUXES COOLING RATES 

WATTS PER SQUARE METER DEGREES / DAY 


L 

UP15 

UP9 

DN15 

DN9 

NT15 

NT9 

CR15 

CR9 

1 

288.2 


1.2 


287.0 


6.81 


2 

288.1 


1.8 


286.2 


6.31 


3 

287.7 


2.8 


28A.9 


. A 1 


A 

287.I 


A.O 


283.1 


-•83 


5 

286.5 


5.7 


280.8 


1.66 


6 

286.2 


7-9 


278.3 


0.68 


7 

286.7 

286.3 

10. A 

8.2 

276.3 

278.1 

0.22 

0.32 

8 

289.2 

289.2 

1 A . 2 

12.o 

275-0 

276.3 

0.1A 

0.20 

9 

291.9 

291.9 

17.8 

16. « 

27A.1 

275.0 

0.36 

0.38 

10 

296.5 

296.5 

25. A 

2A .6 

271.1 

271.9 

1 .86 

1.87 

11 

307.5 

307.5 

56.6 

55.9 

250.9 

251.6 

1.98 

1.98 

12 

32A.8 

32A.8 

101.2 

100.6 

2*'3.6 

22A.2 

1 .87 

1.87 

13 

3A8.5 

3A8.5 

156.8 

156.2 

151.7 

192.3 

1.95 

I.96 

1A 

380.3 

380.3 

229.3 

228.8 

151.0 

151.6 

2.01 

2.01 

15 

A 1 A . 8 

A1A.8 

31A.0 

313.6 

100.8 

101.2 

3.25 

3.26 

16 

A23-5 

A23.5 

3A2.0 

3A1.5 

81.5 

81.9 





TABLE 4 


MIDLATITUDE WINTER PROFILE 


L 

PA 

PU 

TA 

Qt 

03 

» 

• 

1.4 

1.9 

258.1 

0.400E-05 

0.779E-05 

2 

2.7 

3.7 

243.0 

0.400E-05 

0.113E-04 

3 

5.2 

7.2 

227.8 

0.400E-05 

0.116E-04 

k 

10.0 

13.9 

220.6 

0.400E-05 

0.108E-04 

5 

19-3 

26.8 

218.4 

0.400E-05 

0.921E-05 

6 

37-3 

51.8 

215.3 

0.400E-0> 

0.691E-05 


72.0 

100.0 

215.7 

0.400E-05 

0.355E-05 

6 

125.8 

152.9 

217.5 

0.446E-05 

0.165E-05 

9 

187.7 

224.1 

218.7 

0.902E-05 

0.87 IE -06 

10 

269.2 

516.1 

221.6 

0.223E-04 

0.347E-06 

n 

373.6 

433.3 

234.7 

0.986E-04 

0.150E-06 

12 

504.6 

578.3 

247.5 

0.410E-03 

0.848E-07 

13 

665.4 

755 3 

259.8 

0. 1 10E-02 

0.550E-07 

\k 

860. t 

968.0 

267.6 

0.200E-02 

0.467E-07 

15 

992.9 

1018.0 

271.5 

0.257E-02 

0.462E-07 


SURFACE TEMPERATURE - 272.2 


FLUXES COOLING RATES 

WATTS PER SQUARE METER DEGREES / DAY 


L 

UP15 

UP9 

0N15 

DN9 

NT15 

NT9 

CR15 

CR9 

1 

236.4 


1.0 


235.4 


5.10 


2 

236.2 


1.3 


234.9 


3.99 


3 

235.9 


1.9 


234.0 


2.20 


4 

235.6 


2.5 


233-1 


1.45 


5 

235-5 


3.5 


232.0 


1.19 


6 

235.6 


5-5 


230.1 


0.78 


7 

236.3 

236.3 

8.5 

8.0 

227.8 

228.3 

0.49 

0.53 

8 

238.2 

238.2 

13.2 

13.0 

225.0 

225.3 

0.39 

0.39 

9 

240.3 

240.3 

17.7 

1/ .4 

222.6 

222.8 

0.45 

0.45 

10 

243.1 

243.1 

24.4 

24.1 

218.8 

219.0 

0.67 

O.67 

11 

248.1 

248.1 

36.7 

36.5 

211.5 

211.7 

1.33 

1.33 

12 

258.3 

258.3 

65.3 

65.1 

193.0 

193.2 

1.68 

1.68 

13 

273.9 

273.9 

109.7 

109.5 

164.2 

164.4 

1.47 

1.47 

14 

291.9 

291.9 

158.5 

158.3 

133.4 

133.6 

1.07 

1.07 

15 

306.9 

306.9 

200.5 

200.3 

106.4 

106 6 

1.53 

1.53 

16 

311.5 

3H.5 

214.2 

214.0 

97.3 






TABLE 5 


SUBARCTIC SIMMER PROFILE 


L 

PA 

PU 

TA 

QL 


03 

1 

1 .6 

1.9 

277.5 

0.600E- 

■05 

0.570E-05 

2 

2.7 

3-7 

269.6 

0.600E- 

-05 

0.979F-05 

3 

5-2 

7-2 

255.7 

0.600E 

-05 

0.132E-06 

6 

10.0 

13.9 

263.3 

0.800E- 

-05 

0.119E-06 

5 

19-3 

28.8 

238.6 

0.600E 

-05 

0.811E-05 

6 

37-3 

51.8 

225.1 

0.600E- 

-05 

0.519E-05 

7 

72.0 

100.0 

225.1 

0.600E- 

-05 

0.368E-05 

8 

125.6 

152.5 

225.0 

0.801E- 

-05 

0.163E-05 

9 

186.9 

223.0 

225.1 

0 - 1 1 1 E- 

■08 

0.787E-06 

10 

267.6 

316.1 

225-0 

0.622E- 

-06 

0.316E-06 

n 

371.1 

830.2 

260.7 

0.292E 

-03 

0.165E-06 

12 

500.8 

573-9 

256.0 

0.102E- 

■02 

0.997E-07 

13 

680.2 

789.3 

268.7 

0.25U- 

-02 

0.689E-07 

16 

853.1 

960.0 

279.6 

0.689E- 

-02 

0.515E-07 

15 

988.9 

1010.0 

286.0 

0.697E- 

-02 

0.618E-07 


SURFACE TEMPERATURE - 287.0 


FLUXES COOLING RATES 

WATTS PER SQUARE METER DEGREES / DAY 


L 

UP15 

UP9 

DN15 

DN9 

NT 15 

NT9 

CR15 

CR9 

1 

270.3 


1.3 


269.0 


7.18 


2 

270.1 


1-9 


268.2 


6.93 


3 

269.7 


3.0 


266.7 


6.85 


6 

269.0 


6.3 


266.7 


2.91 


5 

268.3 


5-9 


262.6 


1.83 


6 

267.8 


8.2 


259.6 


0.97 


7 

268.0 

268.0 

11.3 

10.0 

256.8 

258.0 

0.55 

0.65 

8 

270.2 

270.2 

16.6 

15.9 

253.6 

256.3 

0.39 

0.61 

9 

272.3 

272.3 

21.1 

20.6 

251.2 

251.7 

0,51 

0.53 

10 

275.0 

275.0 

28.1 

27.7 

266.9 

267.6 

1.01 

1.02 

11 

280.6 

280.6 

66.7 

66.3 

236.0 

236.6 

1.93 

1.93 

12 

296.3 

296.3 

86.8 

86.6 

209.5 

209.8 

1.90 

1.90 

13 

318.8 

318.8 

161 .6 

161. 3 

177.2 

177.5 

1.70 

1.70 

16 

366.0 

366.0 

206.2 

203.9 

161 .8 

162.1 

1.65 

1.65 

15 

376.8 

376.8 

27».l 

270.8 

105.7 

106.0 

2.26 

2.27 

16 

386.9 

386.9 

292.7 

292.6 

92.3 

92.5 





TABLE 6 


SUBARCTIC WINTER PROFILE 


L 

PA 

PU 

TA 

QL 

03 

1 

1 .k 

1-9 

265-9 

0.600E-05 

0.932E-05 

2 

2.7 

3.7 

236.2 

0.600E-05 

0.126E-06 

3 

5-2 

7.2 

222.7 

0.600E-05 

0.121E-06 

6 

10.0 

13.9 

216.5 

0.600E-05 

0.920E-05 

5 

19.3 

26.8 

216. 3 

0.600E-05 

0.87 IE-05 

6 

37.3 

51.8 

213-1 

0.600E-05 

0.780E-05 

7 

72.0 

100.0 

215.7 

0.600E-05 

0.528E-05 

8 

125.7 

152.6 

217.3 

0.509E-05 

0.269E-05 

9 

18;.2 

223.6 

217.3 

0.989E-05 

0.129E-05 

10 

268.2 

316.9 

217-3 

0.169E-06 

0.628E-06 

11 

372.1 

631.6 

225.8 

0.663E-06 

0.131E-06 

12 

502.2 

575.5 

239.6 

0.272E-03 

0.669E-07 

13 

662.1 

751.6 

251.7 

0.762E-03 

0.676E-07 

It 

855.7 

963.0 

258.3 

0.997E-03 

0.355E-07 

15 

987.9 

1013.0 

257.5 

O.898E-03 

0.307E-07 


SURFACE TEMPERATURE « 257-1 


FLUXES COOLING RATES 

WATTS PER SQUARE METER OEGREES / DAY 


L 

UP15 

UP9 

DN15 

DN9 

NT15 

NT 9 

CR15 

CR9 

1 

202.2 


0.8 


201.6 


3.90 


2 

202.1 


1.1 


201.0 


3.22 


3 

201.9 


1.6 


200.3 


1.98 


6 

201.6 


2.1 


.199.5 


1.35 


5 

201.5 


3-1 


198.6 


1.15 


6 

201.6 


6.9 


196.7 


0.86 


7 

202.2 

202.6 

8.1 

8.2 

196.1 

196.2 

••63 

0.65 

8 

203.7 

203.7 

13.1 

13.2 

190.5 

190.5 

0.69 

0.68 

9 

205.1 

205-1 

17.6 

17.6 

187.5 

187.5 

0.51 

0.51 

10 

207.0 

207.0 

23.9 

23.8 

183.2 

183.2 

0.52 

0.52 

11 

210.2 

210.2 

32.7 

32.7 

177.5 

177.6 

0.82 

0.82 

12 

217.3 

217.3 

51.1 

51.0 

166.2 

166.2 

1.57 

1-57 

13 

230.5 

230.5 

91.2 

91.1 

139.3 

139-6 

1.36 

1.36 

16 

263-9 

263-9 

133.1 

133.0 

110.9 

m.o 

O.78 

0.78 

15 

266.5 

266.5 

155*1 

155.1 

91.6 

91.6 

1.06 

1.06 

16 

268.0 

268.0 

162.8 

162.7 

85.2 

85.2 
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Levels in the Radiation Model 
FIG. 1 
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Schematic showing partial cloudiness in one layer 

FIG. 2 
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Schematic showing maximum overlap 
FIG. 3 
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APPENDIX 


CODE LISTING 

NOTE: Several loops diagnosed by the compiler as being non-vectorizable can in 

fact be vectorized by using vector functions or slight restructuring. Since 
these are usually system-specific, the user is expected to make the necessary 
changes. A fully vectorized half-precision version of the code designed for 
the CYBER-205 that is currently being used in the GCM is available from the 
authors upon request. 
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0001 


0002 

0003 

0001 * 

0005 

0006 

0007 

0008 

0009 

001C 

0011 
0012 
0013 
00 u» 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
0021 * 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 
0033 
0031* 


SUBROUTINE RAD IR (LM,PL,PL20,TL,TG,QL,03L,CLD,CLU, 
1 ATLX , RS , ULVTTOP , DIWBOT) 

C 


c********** 

r 

INPUT PARAMETERS *********** 

c 

LM 

..NUMBER OF LAYERS 

c 

PL 

..LAYER MIDDLE PRESSURE IN MB 

c 

PL20. . 

..LAYER BOTTOM PRESSURE IN MB 

c 

TL.... 

..LAYER MIDDLE TEMPERATURE IN KELVIN 

c 

TG 

..SURFACE TEMPERATURE IN KELVIN 

c 

QL.... 

..WATER VAPOR MIXING RATIO IN G/G 

c 

03L... 

..OZONE MIXING RATIO IN G/G 

c 

CLD... 

..RANDOM OVERLAP CLOUD FRACTION 

c 

CLU... 

..MAXIMUM OVERLAP CLOUD FRACTION 


C 

C********** OUTPUT PARAMETERS ********** 

C 

C ATLX LAYER COOLING RATE IN CELSIUS/DAY 

C RS NET UPWARD FLUX AT SURFACE IN W/M**2 

C ULWTOP... UPWARD FLUX AT TOP OF ATMOSPHERE IN W/M**2 
C DLWBOT... DOWNWARD FLUX AT SURFACE IN W/M**2 
C 

DIMENSION PL (72,20) ,PL20 (72,20) ,TL (72,20) ,TG(72) 

DIMENSION QL (72,20) ,03L (72,20) ,CLD (72,20) ,CLU (72,20) 

D I MENS I ON ATLX (72 , 20) , RS (72) , ULWTOP (72) , DLWBOT (72) 

DIMENSION FE (72,20) .COOLR (72,20) 

COMMON/FLUX/FD (72,20) ,FU(72,20) 

COMMON/AMHN/TA (72,20) ,PA (72,20) ,WA (72,20) ,0A (72,20) 
l.NDATA.TS (72) 

COMMON/CSC/CSS (72 , 20) , CCU (72 , 20) , CH (72 , 20) , CC (72 , 20) 

1 , CS (72 , 20,20) 

COMMON/PR tSS/PU (72 , 20) , DP (72 , 20) 

COMMON/ I N I T/NP , NPP 1 , NPM1 
C 

LOGICAL FIRST/. TRUE./ 

IF (.NOT.F IRST) GO TO 10 
F IRST-. FALSE. 

CALL INDATA 
C 

10 CONTINUE 
NDATA-LM 

FAC-3- 6*21*. /10030.*. 98 
DO 20 1-1,72 
PU (I , l)-0. 

PU (I , 2) -1 . 

cssd , 1 ) -0 . 
ccud , 1 ) «o . 

TS(I)-TG(I) 

20 CONTINUE 

DO 30 J-l.NDATA 
DO 30 1-1,72 
PA (I , J) -PL (I , J) 

PU (I , J+2) -PL20 ( I , J) 

TA ( I , J) -TL (I , J) .'R .CEDING PAGE BLANK NOT FILMED" 

WA (I , J) -QL (l , J) 

OA ( I , J) -03L ( I , J) 

CSS(I ,J+1)-CLD (I ,J) 

CCU(I ,J+1)-CLU(I ,J) 

CSS(I , J+l) -AMAX1 (AMIN1 (CSS ( I ,J+1) , 1 .0) ,0.) 



0035 


CCU ( 1 , J+0-AMAX1 (AMIN1 (CCU(l ,J+1) ,1.0) ,0.) 

0036 

30 

r 

CONTINUE 

0037 

L 

r 

CALL SETUP 

0038 

l 

00 40 IP-l.NP 

0039 


DO 40 IL-1,72 

0040 


FE(IL.IP) -FD(IL, IP) -FU (l L, IP) 

0041 

40 

CONTINUE 

0042 


DO 50 1-1,72 

0043 


ULWTOP (l)-FU(l ,1) *0.001 

0044 


DLWBOT (l)-FD(l,NP) *0.001 

0045 


RS(I)— FE(I,NP) *0.001 

0046 

50 

CONTINUE 

0047 


DO 60 IP-2.NP 

0048 


DO 60 IL-1,72 

0049 


COOLR ( 1 L , 1 P) - (FE ( 1 L , 1 P) -FE ( 1 L , 1 P- 1 ) ) *F AC/DP ( 1 L , 1 P) 

0050 

60 

CONTINUE 

0051 


DO 70 J-l.NDATA 

0052 


DO 70 1-1.72 

0053 


ATLX(I , J) -COOLR (1 ,J+2) 

0054 

70 

CONTINUE 

0055 


RETURN 

0056 


END 


UNCOLLAPSIBLE LOOPS REASON FOR NON-VECTOR I ZAT I ON 

LINE 0026 LINE 0035 STATEMENT CONTAINS NONVECTOR I ZABLE F UNCTION 

LINE 0025 LINE 0026 NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 

UNE 0018 LINE 0019 LHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 

NUMBER OF LOOPS IN ROUTINE - 0010 

NUMBER OF VECTOR I ZABLE LOOPS = 0007 

NUMBER OF STACKL I BEO LOOPS * 0000 


0001 SUBROUTINE INDATA 
C 

C ********* READ IN DATA *********** 

C 

0002 C0MM0N/H20/GL (30,20) ,COEFF (30,20,2) 

0003 COMMON/ I N0AT/NUP1 ,NW 

0004 NUP1-20 

0005 NW-22 

0006 READ (5. 10) DUMMY 

0007 10 FORMAT (A 1 , / , A 1 , /) 

0008 READ (5, 10) DUMMY 

0009 READ (5,20) ((GL (IW, IU) , IW-l.NW) , IU-1.NUP1) 

0010 READ (5, 10) DUMMY 

0011 READ (5,20) ( ( (COEFF (I W, I U, I ) , I W-l ,NW) , I U-l ,NUP1) , I -1 , 2) 

0012 20 FORMAT (1 1F7.L) 

0013 RETURN 

0014 END 
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UNCOLLAPSIBLE LOOPS 


REASON FOR NON-VECTORIZATION 


LINE 0011 

LINE 0011 
LINE 0009 

NUMBER OF 

NUMBER OF 

NUMBER OF 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 

0012 

0013 

0014 

0015 
0016 

0017 

0018 

0019 

0020 


0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 


LINE 0011 LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0011 STATEMENT IS NOT AN ASSIGNMENT STATEMENT 

LINE 0009 STATEMENT IS NOT AN ASSIGNMENT STATEMENT 


LOOPS IN ROUTINE - 0003 
VECTOR I ZABLE LOOPS - 0000 
STACKLIBED LOOPS - 0000 


SUBROUTINE SETUP 
C 
C 

DIMENSION PRE (2) ,PR(2) 

DIMENSION X2 (72 , 20) , X4 (72 , 20) . SZ 1 (72 , 20) , SZ2 (72 , 20) , XX (72 , 20) 
DIMENSION RV (72 , 20) , RW (72 , 20) , PM (72 , 20) , PMV (72 , 2Cl , PMW (72 , 20) 
DIMENSION DPI (72,20) ,PAI (72,20) ,TAI (72,20) ,RAWI (72,20) 

* , OZAI ( 72 , 20 ) 

COMMON/PRESS/PU (72 , 20) , DP (72 , 20) 

COMMON/TEMP/TUI (72,20) 

COMMON/ I N I T/NP , NPP 1 , NPM1 

C0MM0N/SH20/UBAR (72,20) ,VBAR (72,20) ,WBAR (72,20) .UBARM (72,20) , 
1 VB ARM (72, 20) , WB ARM (72, 20) 

C0MMON/SC0203/PSCALV (72 , 20) , PSCALW (72 , 20) , U (72 , 20) , V (72 , 20) , 

1 PSCAMV (72 , 20) , PSCAMW (72 , 20) , UM (72 , 20) , VM (72 , 20) 

COMMON/AMHN/TA (72 , 20) , PA (72 , 20) , WA (72 , 20) , OA (72 , 20) , NDATA 
1 . TS (72) 

DATA PRE/275 • . 550 . / , PR/30 . , 300 . / , PP/ 1 . / 

NP-NDATA+2 
NPP1-NP+1 
NPM1-NP-1 
DO 10 IP-2.NP 
DO 10 .1 L-l ,72 

DP ( I L , I P) -PU ( I L , IP) — PU ( I L , I P— 1 ) 

DPI ( I L , I P) "DP ( I L , I P) * 1 . 02 
10 CONTINUE 
C 

C ****** TEMPERATURE AND HUMIDITY INTERPOLATIONS ************* 

C 

DO 20 J-l, NDATA 
DO 20 1-1,72 

WA (I , J) -AMAX1 (WA ( I , J) ,0. IE-12) 

OA ( I , J) -AMAX1 (OA (I , J) ,0. IE-12) 

20 CONTINUE 

DO 30 J-3.NP 
DO 30 1-1,72 
PA I (I , J) -PA (I , J-2) 

TAI (I ,J) -TA (I , J-2) 

RAW I (I , J)-WA (I , J-2) *DPI (I.J) 

OZAI (I , J)-OA (I , J-2) *DPI (I.J) 

30 CONTINUE 

DO 40 J-3.NPMI 
DO 40 1-1,72 

X2 (I , J) -ALOG (PU ( I , J) /PA ( I , J-2) ) 



0036 

0037 

0038 
0039 
0060 
006) 
0062 

0063 

0066 

0065 

0066 

0067 

0068 
0069 
0050 

0051 

0052 
0053 
0056 


0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 
0063 
0066 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 


0073 

0076 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 
0083 


X6 (I ,J) -AL0G (PA ( I , J- 1) /PA ( I , J-2) ) 

TUI (1 , J)-TA (I » J-2)+X2 (I , J) /X6 (I , J) *(TA (I , J-l) -TA (I , J-2)) 

60 CONTINUE 
10 50 (•1,72 

X2 (I ,NP) -ALOG (PU (I ,NP) /PU (I ,NPM1)) 

X6 (I ,NP) -ALOG (PA (I .NDATA) /PU (I ,NPM1) ) 

TUI (l.NP)-TUI (I ,NPM1)+X2 (I ,NP) /X6(l ,NP)*(TA(I ( NDATA) 

1 -TUI (I.NPM1)) 

TUI (I ,NPP1)-TS (I) 

OPI (l,l)-0. 

PA I (l f l)-0. 

PA I (l,2)-0.5*PU(l,2) 

50 CONTINUE 
DO 60 J-l ,2 
DO 60 1-1.72 
TAI (I , J) -TA (1,1) 

TUI (I.J)-TA(I.I) 

RAW I (I , J)-WA (1 . 1) *DPI (I.J) 

OZAI (I.J) -OA (1,1) *DP I (I.J) 

60 CONTINUE 
C 

C********** COMPUTE SCALED WATER VAPOR AMOUNTS********* 

C 

DO 80 IL-1,72 
UBAR ( I L , 1 ) -0 . 

VBAR(IL, D-0.0 
WBAROL.D-0.0 
80 CONTINUE 

DO 90 IP-2.NP 
DO 90 IL-1,72 

SZ1 (IL, IP) -EXP (0.005* (TAI (IL, IP) -225.0)) 

SZ2 (I L, IP) -EXP (0.016* (TAI (I L, IP) -256.0)) 

VBAR ( I L , I P) -VBAR ( I L , IP-1)+RAWI ( I L . IP) *SZ 1 (IL.IP)* 

1 (PAI (IL.IP) /PRE(l)) 

WBAR (IL.IP) -WBAR (I L, IP-D+RAWI (I L, IP) *SZ2 (I L, IP) * 

1 (PAI (IL.IP) /PRE (2)) 

XX (IL.IP) -(PAI (IL.IP) /630 . ) *RAW I (IL.IP) *RAW I (IL.IP) /DPI (IL.IP) 
XX (I L, IP) -XX (IL.IP) *EXP (1800./TAI (IL, IP) -6.0811) 

UBAR(IL, IP) -UBAR(IL. IP-D+XX (IL.IP) 

UBARM (I L, IP) -0.5* (UBAR (I L, I P)+UBAR ( I L , I P-1) ) 

VBARM ( I L . I P) -0 . 5* (VBAR (IL.IP) +VBAR ( I L . I P- 1 ) ) 

WBARM (IL.IP) -0 . 5* (WBAR (IL.IP) +WBAR ( I L , I P- 1 ) ) 

90 CONTINUE 
C 

C********* COMPUTE SCALED n2 AMOUNTS ********** 

C 

DO 100 IL-1,72 
PSCALV (I L . 1) -0. 

PSCALW(I L, l)-0. 

100 CONTINUE 

DO 110 IP-2, NP 
DO 110 IL-1,72 

RV(IL, I P) -EXP (O.OO89* (TAI (IL.IP) -260. 0 ) ) 

RW ( I L, IP) -EXP (0.025* (TAI (IL, I P) —260 .C) ) 

PM(IL, IP) -AMAX1 (PAI (IL.IP) ,PP) 

PMV (IL.IP) -(PM (IL.IP) /PR (1))**0.85*RV (IL.IP) 

PMW (IL.IP)- (PM (IL.IP) /PR (2) ) **0.5*RW(IL, IP) 

C 

C********** 0.26 IS FOR C02-330 PPMV *********** 

C 



0084 

0085 

0086 

0087 

0088 


0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 


PSCALV (I L • IP)-DP(IL, IP)*PMV(IL, IP)*0.26+PSCALV(IL, IP-1) 
PSCALW(IL, IP)-DP(IL, IP)*PMW(IL, IP)*0.26+PSCALW(IL, IP-1) 
PSCAMV ( I L , I P) -0 . 5* (PSCALV ( I L , I P) +PSCALV ( I L , I P- 1 ) ) 
PSCAMW ( I L , I P) -0 . 5* (PSCALW ( I L , I P) +PSCALW ( I L , 1 P- 1 ) ) 

110 CONTINUE 
C 

C ********** COMPUTE SCALED OZONE AMOUNTS ******** 

C 

DO 120 IL-1,72 
U(IL,l)-0. 

V(IL, l)-0. 

120 CONTINUE 

DO 130 IP-2, NP 
DO 130 IL-1,72 

U ( I L , I P) -PA I ( I L , I P) *OZA I ( I L , I P) +U ( I L , I P- 1 ) 

V ( I L , I P) -OZA I ( I L , I P) +V ( I L , I P - 1 ) 

UM ( I , I P) -0 . 5* (U ( I L , I P) +U ( I L , I P- 1 ) ) 
VM(IL,IP)-0.5*(V(IL,IP)+V(IL,IP-1)) 

130 CONTINUE 
C 

CALL IRFLUX 
C 

RETURN 

END 


UNCOLLAPSIBLE LOOPS REASON FOR NON-VECTORIZATION 


RHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 
RHS ARRAY h.;i POSSIBLY RECURSIVE PROPERTIES 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 
RHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 
PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR I ZAT I ON OF LOOP 

LHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 
STATEMENT CONTAINS NONVECTOR I ZABLE FUNCTION 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 

NUMBER OF LOOPS IN ROUTINE - 0020 
f. JMBER OF VECTOR I ZABLE LOOPS - 0010 
NUMBER OF STACKLIBED LOOPS - 0000 


LINE 0094 
LINE 0093 
LINE 0078 
LINE 0077 
LINE 0061 
LI»JE 0060 
LINE 0048 

LINE 0039 
LINE 0022 
LINE 0021 


LINE 0096 
LINE 0094 
LINE 0085 
LINE 0078 
LINE 0068 
LINE 0061 
LINE 0053 

LINE 0045 
LINE 0024 
LINE 0022 


0001 SUBROUTINE IRFLUX 
C 

C 

0002 DIMENSION SV(7?., 20) ,SW(72, 20) , TU1 1 (72,20; ,TUI2 (72.20) 

0003 DIMENSION BAI (72.20) ,BWI (72,20) 

0004 D I MENS ION X 1 (72 , 20) , X2 (72 . 20) , X3 (72 . 20) , X4 (72 , 20) , X5 (72 , 20) , 

1 X6 (72 , 20) , X7 (72 , 20) , X8 (72 , 20) 

0005 DIMENSION BH20VC25) .BH20W(25) .BLKC02(25) ,BLKWIN(25) 

0006 DIMENSION FH (72,20) ,DH (72,20) ,F1 (72,20) ,IT(72,20) 

0007 DIMENSION SH (72,20) ,SG (72,20,20) ,SHT (72,20) , SHU (72,20) ,SUM(72) 

0008 DIMENSION A ( 0001 0002 0003 0004 0005 0006 7 0008 2,20) ,B (72,20) ,C (72,20) ,WV (72, *0) ,YV(72,20) , 



0009 

0010 

0011 

0012 

0013 

0014 

001 $ 

0016 


0017 


0018 


0019 


0020 


0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0u40 


1 ZV (72,20) ,WW(72,20) ,YW(72,20) ,ZW(72,20) 

D I MENS I ON D (72 . 20) , E (72 ,20) f F (72 1 20) , TAUF (72 , 20) , TXUF (72 , 20) 
C0MM0N/SH20/UBAR (72 , 20) , VBAR (72 , 20) , WBAR (72 , 20) , UBARM (72 , 20) , 

1 VBARM (72 , 20) , WBARM (72 , 20) 

COMMON/SC02O3/PSCALV (72,20) ,P$CALW(72,20) ,U (72,20) ,V(72,20) , 

1 PSCAMV (72 , 20) , PSCAMW (72 , 20) , UM (72 , 20) , VM (72 , 20) 

COMMON/ I N I T/NP , NPP 1 , NPM 1 
COMMON/FLUX/FD (72,20) ,FU (72,20) 

COMMON/CSC/CSS (72 , 20) , CCU (72 , 20) , CH (72 , 20) , CC (72 , 20) 

1, CS (72,20,20) 

COMMON/1 EHP/TU I (72,20) 

DATA TEMPI/ 190./, DT/5./,NT/25/ 

C 

C******* H20 PLANCK FUNCTION ********** 

C 

DATA BH20V/2266 1 . , 24086 . , 2;>575 . . 27 1 35 . . 28775 . . 30506 . , 32339 . . 

* 34286 . , 3636 1 . , 38578 . , 40954 . , 43505 . , 46248 . , 49203 • , 52388 . , 

* 55824 . , 59532 . , 63533 • , 67849 . , 72502 . , 775 1 6 . , 829 1 3 • . 887 1 7 . . 

* 94952., 101640./ 

DATA BH20W/30708 . , 34405 . , 384 1 7 . , 42763 . , 4746 1 . , 52529 • . 57985 . . 

* 63850.. 70141., 76880., 84088., 91784., 99990., IO8726., 118016., 

* 127881 . , 138344. , 149429. , 161160. ,17356 1 . , 186659. ,200478. , 

* 2 15046., 230390., 246539-/ 

C 

C******* WINDOW PLANCK FUNCTION ********** 

C 

DATA BLK W I N/ 1 944 . , 2377 . . 2877 • . 3450 . , 4 1 02 . , 4838 . , 5664 . , 6585 . . 

* 7606 . , 8733 • . 9969 . , 1 1 320 . , 1 1788 . , 1 4380 . , 1 6097 • . 

* 1 7944 ..19923.. 22038 . , 24292 . , 26685 . , 2922 1 . , 3 1 902 . , 

* 34729.. 37703., 40825-/ 

C 

C******* C02 PLANCK FUNCTION ********* 

C 

DATA BLKC02/ 1 8606 .,21145., 23884 . , 26826 . . 29973 • . 33325- » 36883 . . 

* 40647 .,44617., 48792 . , 1 70 . ,57750 . , 62530 . , 67509 • . 

* 72683 . , 78050 . , 83609 . , 89354 • , 95285 ..101397., 107688 . , 

* 1 14155. . 120794. , 127601 . , 134574./ 

DO 10 IP-1, NP 

DO 10 IL-1,72 
SHT(IL, IP) -0. 

SHU (I L, I P) -0. 

10 CONTINUE 

DO 15 IP-1, NP 
DO 15 IX-1.NP 
DO 15 IL-1,72 
SG (I L, IX. IP) -0. 

15 CONTINUE 

DO 20 IP-1, NPP1 
DO 20 IL-i,72 

FH(IL, IP) -(TUI (IL, IP) -TEMP 1)/DT+ 1.5 
IT (I L, IP) -FH (IL, IP) 

I T ( I L , I P) -MAXO (M I NO ( I T ( I L , IP) ,NT-1) ,1) 

F 1 (I L, I P) -FLOAT (IT (I L, I P) -1) 

DH (I L, I P) -TUI ( I L , I P) - (TEMP 1+F 1 ( I L , I P) *PT) 

SV(IL.IP) -BH20V ( I T ( I L , I P) ) + (BH20V ( I T ( I L , IP) + 1) -BH20V (IT (I L, IP) 

1 ))*OH(IL,IP)/DT 

SW(IL,IP) -BH20W(IT (I L, I P) ) + (BH20W (IT (I L, IP)+1) -BH20W(IT (I L, IP) 

1 ))*DH(IL,IP)/DT 

BA I (IL, IP) -BLKC02 (IT (IL, IP) ) + (BLKC02 (IT(IL,IP) + 1) -BLKC02 (IT 
1 ( I L , I P) ) ) *DH ( I L , I P) /DT 



20 


0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 


0049 

0050 

0051 

0052 

0053 

0054 

0055 


0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 
0066 

0067 

0068 
0069 
0070 
0071 


0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 


30 


BWI (It, IP) -BLKWIN (IT (I L, IP) ) + (BLKWIN (IT (I L, IP)+1) -BLKWIN (IT 
1 ( I L , I P) ) ) *DH ( I L , I P) /DT 

CONTINUE 
DO 30 I P-1, NPP 1 
DO 30 11-1.72 

SH ( I L , I P) -SV ( I L , I P) +SW ( 1 1 , 1 P) +BA I (IL,IP)+BWI (IL.IP) 

TU 1 1 ( I L , I P) -TU I (IL.IP) -250. 

TUI2(IL, IP)-TUI 1 (IL.IP) *TUI1 (IL.IP) 

CONTINUE 


C ***************** H20 **************** 

C ***************** **************** 


c 


40 


C 

C 

C 


45 


50 


C 

c 

c 


DO 40 IP-2.NP 
DO 40 5 L-l ,72 

A (IL, IP-1)-ABS (UBAROL, IP)) 

B (IL, IP-1)-ABS (VBAR (IL.IP)) 

C (IL, IP-1)»ABS(WBAR(IL, IP)) 

CONTINUE 

JP-NPH1 

CALL TABLE (JP,A,B,C,WV,YV,ZV,WW,YW,ZW) 

CALL GFUNC (JP.A.B.C.WV, YV.ZV.WW.YW.ZW) 

DO 45 IP-2.NP 
DO 45 IL-1./2 

XI (IL.IP) -WV ( I L , I P- 1 ) *SV ( I L , 1 ) 

X2 ( ! L , I P) - 1 .+YV ( I L , I P- 1 ) *TU 1 1 (I L, 1)+ZV(I L, I P-1) *TUI 2 (I L, 1) 
X3(IL, IP)-WW(IL, IP-1)*SW(IL, 1) 

X4(IL,IP)-1.+YW(IL, I P-1) *TUI 1 (I L, l)+ZW(l L, I P~l) *TUI 2 (I L, 1) 
SHT(IL, IP) -XI (IL.IP) *X2 (I L , IP) +X3 (I L , IP) *X4 (IL.IP) 

CONTINUE 

DO 50 IP-l.NPMl 

DO 50 IL-1.72 

A (I L . I P) -ABS (UBAR (I L . NP) -UBAR (IL.IP)) 

"'(IL.IP) -ABS (VBAR(IL.NP) -VBAR (IL.IP)) 

L (I L , 1 P) -ABS (WEAR (I L , NP) -WBAR (IL.IP)) 

CONTINUE 

JP-NPMl 

CALL TABLE (JP.A.B,'" WV, YV.ZV.WW.YW.ZW) 

CALL GFUNC (JP,A.B,C,WV,YV,2V.WW. YW.ZW) 


DO 55 IP-l.NPMl 
DO 55 IL-1.72 

XI (IL, IP)-WV(IL. IP) *SV (I L.NPP1) 

X2 ( I L , I P) - 1 .+YV ( I L , IP) *TUI 1 ( I L . NPP 1 ) +ZV (IL.IP) *TU 1 2 ( I L , NPP 1 ) 
X3 (I L. IF)-WW(IL, IP) *SW (I L.NPP1) 

X4(IL, IP)-1.+YW(IL, I P) *TU 1 1 (I L,NPP1)+ZW (I L, IP) *TUI2 (IL.NPP1) 
X5 ( I L , IP) -WV (IL.IP) *SV (I L.NP) 

X6 (I L. I P) -1 ,+YV (I l , I P) *TUI I ( I L , NP) +ZV (IL.IP) *TU 1 2(1 L.NP) 

X7 (I L . I P) -WW (IL.IP) *SW ( I L , NP) 

X8(IL, IP)-1.+YW(IL, IP) *TUI 1 ( I L , NP) +ZW (IL.IP) *TU 1 2 ( I L , NP) 

SHU (I L, I P) -XI (IL.IP) *X2 (I L, IP) +X3 (IL.IP) *X4 (I L, I P) - 
1 X5 ( I L , I P) *X6 (IL.IP) -X7 (IL.IP) *X8 (IL.IP) 

55 ‘CONTINUE 


C 


C*********DOWNWARD FLUX SG (I X , I P) ********** ** 


C 



60 


0085 

0086 

0087 

0088 
0089 
0090 

0091 

0092 


0093 

0091 * 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 


0107 

0108 

0109 

0110 
0111 
0112 

0113 

0114 

0115 


Cl 16 

01 17 
01 18 
0119 
0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 
0129 


C 

C 

C 


00 75 IP-2, HP 
00 60 IX-2.IP 
00 60 I 1-1,72 

A( I L . I X- 1 ) -AOS (UBAR ( I L , I P) -UBAIM ( I L . I X) ) 
B (II, I X— 1) — ABS (VBAR ( I L , I P) - VBARM ( I L , I X) ) 
C (I L, IX- 1) -ABS (UBAR (I L, I P) -WBARft (I l, I X) ) 
CONTINUE 
JP-IP-1 

CALL TABLE (JP, A,B,C.WV,YV,ZV,WW, YV,ZW) 
CALL GFUNC (JP, A,B,C,WV, YV,ZV,WW,YW,ZW) 


00 70 IX-2.IP 
DO 70 IL-1,72 

XI (IL, IX-1)-WV(IL. IX-1) *SV(IL, IX-1) 
X2(IL,IX-1)-1.+YV(IL,IX-1)*TUI 1 (IL, IX-1)+ZV(IL, IX-1) 

1 *TUI2(IL, IX-1) 

X3 (I L. IX-1) -WW(IL, IX-1) *SW(IL. IX-1) 
X4(ll,IX-l)«l.+YV(IL.!X-l)*TUI I (IL, IX-1)+2W(IL, IX-1) 

1 *TUI2(IL, IX-1) 

X5 0 L, IX-1) -WV(IL, IX-1) *SV(IL, IX) 

X6 ( I L , IX-1)»1.+YV(IL, IX-1) *TU 1 1 (IL, IX)+ZV(I L, IX-1) *TUI2 (IL, IX) 
X7 (I L, IX-1) -WW(IL, IX-1) *SW(I L, IX) 

X8(IL, IX-1)-1.+YV(IL, IX-1) *TUI 1 (IL.IX)+ZW(IL.IX-1)*TUI2(IL,IX) 
SG ( I L , I X- 1 , 1 P) -X 1 ( I L , IX-1) *X2 (I L, I X-l)+X3 (I L, IX-1) *X4 (IL, IX-1) 
1 -X5(IL, IX- I) *X6 (I L, IX-1) -X7 (I L, IX-1) *X8 (IL, IX-1) 

70 CONTINUE 
75 CONTINUE 
C 

C* ** * ** ***UPWAR0 FLUX SG ( I X , I P) ********* 

C 


DO 95 IP-1.NPH1 
IP1-NP-IP 
DO 80 IX-1.IP1 
DO 80 IL-1,72 

A ( I L , I X) -ABS (UBAR ( I L , I P) -UBAR* (IL, IX+IP)) 

B (IL, IX) -ABS (VBAR ( I L , IP) -VBARM ( I L , I X+l P) ) 

•- ill, IX) -ABS (WEAR ( I L , IP) -WBAR«(IL, IX+IP)) 

80 CONTINUE 
JP-IP1 
C 

C CALL TABLE (JP,A,B,C,WV,YV,ZV,WW,YW,ZW) 

CALL GFUNC (JP,A,B,C,WV,YV,ZV,WW,YW,ZW) 

C 

DO 90 >X-1 IP) 

DO 90 IL-1,72 

XI (IL.IX)-WV(IL, IX)*SV{IL, IX+IP) 

X2 (I L, IX) -1 .+YV (IL, IX) *TUI 1 (I L, I X+l P)+ZV (I L, I X) *TUI 2 (I L, IX+IP) 
X^IL, IX)-WW(IL, IX)*SW(IL, IX+IP) 

X4UL, IX)-1.+YW(IL, IX) *TUI 1 ( I L , IX + IP) +2W ( I L , IX) *TU 1 2 ( I L , IX+IP) 
X5(IL, IX)-WV(IL, IX)*SV(lu, IX+IP-1; 

X6(IL, I X) -1 ,+YV (I L, I X) *TU I 1 (IL, I X+l P- 1)+ZV (I L, IX) 

1 *TUI2(IL, IX+IP-1) 

X7 0 L, I X) -WW (1 L, IX)*SW(IL, I X+l P-1) 

X8(IL, I X) -1 .+YW ( I L, I X) *TU II ( I L , I X+ 1 P- 1 ) +ZW ( I L , I X) 

1 *TU 1 2 ( I L , IX+IP-1) 

SG(IL, IX+IP, IP)-X1 (IL, IX)*X2(IL, IX)+X3(IL, IX)*X4(IL, IX)- 
1 X5 (I L. IX)*X6(IL, I X) -X7 (I L, I X) *X8 (I L, I X) 

90 CONTINUE 

95 CONTINUE 42 



0130 

0131 

0132 

0133 

0134 

0135 
0135 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 
0)46 
0167 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 


0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 
0173 



C 

DO 140 IP-2.NP 
DO 140 IL-1,72 

A ( I L , I P- 1) -ABS (PSCALV ( I L , I P) ) 

B(IL, IP-1)-ABS (PSCALW ( I L , I P) ) 

C(IL, IP-1) -ABS (WBAR(IL.IP)) 

D (IL, IP-1) -ABS (UBAR ( I L , I P) ) 

E (IL, I P-1) -ABS (V(l L, IP) ) 

F (IL, I P-1) -ABS (U(IL,IP)) 

140 CONTINUE 
JP-NPM1 
C 

CALL TRC03 (JP.A,B,C,D,E,F,TAUF ,TXUF) 

C 

DO 145 1P-2.NP 
DO 145 IL-1,72 

SHT ( I L , I P) -SHT ( I L , I P) +BA I (IL,1)*TAUF (I L. IP-1) 

1 +BWI (I L. 1) *TXUF (I L, IP-1) 

145 CONTINUE 

DO 150 IP-1.NPH1 
DO 150 IL-1,72 

A ( I L , I P) -ABS (PSCALV ( I L , NP) -PSCALV ( I L , I P) ) 

B ( I L , I P) -ABS (PSCALW ( I L , NP) -PSCALW ( I L . I P) ) 

C ( I L , I P) -ABS (WBAR ( I L , NP) -WBAR ( I L , I P) ) 

D ( I L , I P) -ABS (UBAR ( I L , NP) -UBAR ( I L , I P) ) 

E (I L, IP) -ABS (V ( I L , NP) -V (I L, IP) ) 

F(IL.IP) -ABS (U ( I L , NP) -U ( I L , I P) ) 

150 CONTINUE 
JP-NPM1 
C 

CALL TRC03(JP,A,B,C,D,E.F,TAUF,TXUF) 

C 

DO 155 IP-1.NPH1 
DO 155 IL-1,72 

SHU (I L, IP) -SHU (I L , IP)-*- (BA I (I L.NPP1) -BA I (I L.NP) ) *TAUF (IL, IP) 
1 + (BWI (IL.NPPl)-BWI (I L.NP) ) *TXUF (IL, IP) 

155 CONTINUE 
C 

C*********DOWNWARD FLUX SG ( I X , I P) ************ 

C 

DO 175 IP-2.NP 
DO 160 IX-2,l n 
DO 160 IL-1,72 

A (I L, I X— 1) —ABS (PSCALV (I L, IP) -PSCAHV (IL, IX)) 

B(IL, IX-D-ABS (PSCALW (I L, IP) -PSCAMW ( I L , I X) ) 

C (I L, IX- 1) -ABS (WBAR (I L, I P) -WBARM(IL, IX)) 

0 ( I L . I X- 1 ) -ABS (UBAR ( I L , I P) -UBARM ( I L , I X) ) 

E (IL, IX-D-ABS (V(IL.IP)-VM(IL, IX)) 

F ( I L , I X- 1 ) -ABS (U ( I L , I P) -Uft ( I L , I X) ) 

160 CONTINUE 
JP-IP-1 
C 

CALL TRC03 (JP, A.B, C,D, E , F ,TAUF ,TXUF) 

C 

DO 170 IX-2.IP 

DO 170 IL-1,72 43 



017 * 


0175 

0176 


0177 

0178 

0179 

0180 

0 1 8 1 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 


0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 


0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 
0209 
0210 
0211 
0212 
0213 


>G(IL. IX-1 . IP)-SG (IL, IX-1 , IP) + (BAI (I L, IX-1) -BA I (IL.IX)) 

1 *TAUF(IL,IX-1) + (BWI (IL.IX-l)-BWI (IL.IX)) 

2 *TXUF(IL,IX-1) 

170 CONTINUE 

175 CONTINUE 
C 

C*********UPWARD FLUX SC (IX. IP) ********* 

C 

00 195 IP-1, NPM1 
IP1-NP-IP 
DO 180 IX— 1 , 1 PI 
00 180 11-1,72 

A ( I L , I X) -ABS (PSCALV ( I L . I P) -PSCAMV ( I L , I X+ 1 P) ) 

B ( I L , I X) -ABS (PSCALW(IL, IP) -PSCAMW(IL, IX+IP)) 

C ( I L . I X) -ABS (WEAR ( I L , I P) -WBAR.i ( I L . I X+ 1 P) ) 

0 ( I L , I X) -ABS (UBAR ( I L , I P) -UBARM ( 1 L . I X+ 1 P) ) 

E (IL.IX) -ABS (V (I L, I P) -VM (I L, IX+IP)) 

F (IL.IX) -ABS (U(IL, IP) -UM (I L. IX+IP) ) 

180 CONTINUE 
JP-IP1 
C 

CALI TRC03(JP.A.B,C,D,E,F,TAUF,TXUF) 

C 

DO 190 IX-1, IP! 

DO 190 IL- 1,72 

SG(IL, IX+IP, IP) -SG (I L, IX+IP, I P) + (BA 1 (IL, IX+IP) -BA I (IL. IX+IP-1)) 

1 *TAUF (IL, I X) + (BW I ( I L , IX+IP) -BW I (IL, IX+IP-1)) 

2 *TXUF (IL.IX) 

190 CONTINUE 

195 CONTINUE 
C 

CALL CLOUDS 
C 

DO 210 IP-1, NP 
DO 210 IL-1,72 
FU(IL, IP) -0. 

FD (I L, I P) -0. 

210 CONTINUE 
C 

C*********** COMPUTE DOWNWARD FLUXES ********** 

C 

DO 230 IP-2.NP 
DO 215 IL-1,72 
SUM (I L) -0. 

215 CONTINUE 
IPM1-IP-1 
DO 220 IX-1, I PM1 
DO 220 IL-1,72 

SUM ( I L) -SUM ( I L) +SG ( I L , IX, IP) *CS (IL, IX, IP) 

220 CONTINUE 

00 225 IL-1,72 

FD (I L, IP) -SU«() L) -SHT (IL. IP) *CH (I L, IP)+SH (I L. IP) 

225 CONTINUE 
230 CONTINUE 
C 

C********** COMPUTE UPWARD FLUXES ************ 

C 


0214 

0215 

0216 


DO 250 IP-l.NPMl 
DO 235 IL“1,72 
SUM (I L) -0. 


44 



0217 

235 

CONTINUE 

0218 


IPP1-IP+1 

0219 


DO 260 IX-IPP1.NP 

0220 


DO 260 IL-1,72 

0221 


SUM(IL)-SUM(IL)+SG(IL,IX,IP) *CS(IL,IX, IP) 

0222 

260 

CONTINUE 

0223 


DO 265 IL-1.72 

0226 


FU(IL,IP) -SUM(IL)+SHU(IL, IP) *CC (IL, IP)+SH (II, IP) 

0225 

265 

CONTINUE 

0226 

250 

CONTINUE 

0227 


DO 260 11-1,72 

0228 


FU(IL.NP) -SH(IL.NPPl) 

0229 

260 

CONTINUE 

0230 


RETURN 

0231 


END 


UNCOLLAPSIBLE LOOPS 

REASON FOR 

NON- VECTOR 1 2AT 1 ON 

LINE 02’ 6 

LINE 

0226 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0219 

LINE 

0221 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0201 

LINE 

0211 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0206 

LINE 

0208 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0177 

LINE 

0190 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0179 

LINE 

0186 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0160 

LINE 

0172 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0161 

LINE 

0168 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 015^ 

LINE 

0158 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0165 

LINE 

0152 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0161 

LINE 

0163 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0107 

LINE 

0117 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0109 

LINE 

0113 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0085 

LINE 

0096 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0086 

LINE 

O 

cn 

0 

0 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0073 

LINE 

0082 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0065 

LINE 

0069 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0057 

LINE 

0062 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0031 

LINE 

006 V 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0032 

LINE 

0035 

STATEMENT CONTAINS NONVECTOF. 1 ZABLE FUNCTION 

LINE 0026 

LINE 

0027 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED JITHIN LOOP 



NUMBER OF LOOPS IN ROUTINE » 0058 
NUMBER OF VECTOR I ZABLE LOOPS - 0037 
NUMBER OF STACKLIBED LOOPS - 0000 


0001 


0002 

0003 

0004 

0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 
0029 
0030 

0031 

0032 


0033 


0034 

0035 

0036 

0037 

0038 

0039 

0040 


SUBROUTINE TABLE (N.A.B,C,WV,YV,ZV,WW,YW,ZW) 

C 

C 

D I MENS I ON KW (72 , 20) , FW (72 , 20) , KU (72 . 20) , FU (72 , 20) 

01 MENS I ON ADEL (72.20) .BDEL (72.20) ,WDEL (72,20) 

DIMENSION A (72,20) ,B (72,20) ,C (72,20) ,WV (72,20) ,YV(72,20) , 

1 ZV (72,20) ,WW(72,20) ,YW (72,20) ,ZW(72,20) 

COMMON/ I NDAT/NUP 1 , NW 
C0MM0N/H20/GL (30,20) .COEFF (30.20,2) 

DATA SV/- 6 . / , SW/-5 • 4/ , DW/O . 3/ , DU/O . 003/ 

NU-NUP1-1 
DO 10 J-l ,N 
DO 10 1-1,72 

B (I , J) -AMAX1 (B (I ,J) .l.E-6) 

B (I ,J) -AL0G10 (B (I , J) ) 

KW (I , J) - (B (I , J) -SV) /DW+1 .5001 
KW (I , J) -M I NO (KW(I,J) , NW- 1 ) 

F W ( I , J) - (B ( I , J) - (SV+F LOAT (KW ( I , J) - 1 ) *DW) ) /DW 

ADEL (I , J)- (COEFF (KW ( I , J)+l ,NUP1 , 1) -COEFF (KW(I , J) ,NUP1 , 1) ) 

1 *FW (I , J) 

BDEL (I ,J)« (COEFF (KW(I , J) + l ,NUP1 ,2) -COEFF (KW(I , J) ,NUP1 ,2) ) 

1 *FW (I , J) 

WDEL (I , J) - (GL (KW ( I , J) + l ,NUP1) -GL (KW (I , J) ,NUP1))*FW(I ,J) 

WV(I,J) -GL (KW ( I , J) ,NUP1) +WDEL ( I , J) 

WV(I ,J)-AMAX1 (WV ( I , J) .0.0) 

YV (I , J)-0. 1* (COEFF (KW(I , J) ,NUP1 , 1)+ADEL (I , J) ) 

ZV (I ,J) -0.001* (COEFF (KW ( I , J) ,NUP1 ,2)+BDEL (I , J) ) 

KU (I , J) -A (I ,J)/DU+1 .5 
KU (I , J) -Mi NO (KU (I , J) ,NU- 1) 

FU ( I , J) - (A (I ,J) -FLOAT (KU (I , J) - 1) *DU) /DU 
C ( I , J) -AMAX 1 CC Cl • J) J.E-6) 

C (I . J) -AL0G10 (C (I , J) ) 

C ( I , J) -AMAX 1 (C ( I , J) , SW) 

KW(I , J)-(C (I , J) -SW) /DW+1 .5001 
KW ( I , J) -M I NO (KW ( I , J) , NW- 1 ) 

FW (I , J) - (C (I , J) - (SW+FLOAT (KW (I , J) -1) *DW) ) /DW 

ADEL (I , J)- (COEFF (KW ( I , J)+l ,KU (I , J) , 1) -COEFF (KW (I , J) , KU (I , J) , 1) ) 

1 *FW(I ,J) + (COEFF (KW ( I , J) ,KU (I , J) + l , 1) -COEFF (KW (I , J) ,KU (I , J) , 1) ) 

2 *FU (I , J) 

BDE L ( I , J) - (COE F F (KW ( I , J) +1 , KU ( I , J) , 2) -COEF F (KW ( I , J) , KU ( I , J) , 2) ) 

1 *FW(I ,J) + (COEFF (KW ( I ,J) ,KU (I , J) + l ,2) -COEFF (KW(I , J) ,KU(I,J) ,2)) 

2 *FU (I , J) 

WDEL (I , J) - (GL (KW (I , J)+l , KU (I , J) ) -GL (KW (I , J) , KU ( I , J) ) ) *FW (I , J) + 

1 (GL (KW (I , J) , KU (I , J) + 1) -GL (KW ( I , J) ,KU (I , J) ) ) *FU (I , J) 

WW(I ,J)-GL(KW(I ,J) , KU (I , J) ) +WDEL ( I , J) 

WW(I ,J)-AMAX1 (WW ( I ,J) ,0.0) 

YW (I ,J) -0.1* (COEFF (KW (I ,J) .KU (I ,J) , D+ADEL (I , J)) 

ZW (I, J) -0.001* (COEFF (KW (I , J) ,KU(l,J) ,2)+BDEL (I , J) ) 

10 CONTINUE 

RETURN 46 



0041 


END 


UNCOLLAPSIBLE LOOPS REASON FOR NON-VECTOR I ZAT I ON 

LINE 0010 LINE 0038 RHS ARRAY MUST BE REAL, INTEGER, LOGICAL, 

HALF PRECISION. OR COMPLEX 

LINE 0009 LINE 0010 NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 

NUMBER OF LOOPS IN ROUTINE - 0002 
NUMBER OF VECTOR I ZABLE LOOPS - 0000 
NUMBER OF STACKLIBED LOOPS - 0000 
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0006 

0007 

0008 

0009 
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0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 
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0021 
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0024 

0025 

0026 

0027 

0028 
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0030 

0031 


SUBROUTINE GFUNC (N, A,B,C,WV,YV,ZV,WW, YW.ZW) 

C 

C 

DIMENSION A (72,20) ,B(72,20) ,C(72,20) ,WV (72,20) ,YV (72,20) , 

1 ZV (72,20) ,WW(72,20) ,YW(72,20) ,ZW (72,20) 

DIMENSION A1 (72,20) ,A2 (72,20) , A3 (72,20) ,B1 (72,20) ,B2 (72,20) , 

1 B3 (72.20) ,B4 (72,20) . Cl (72,20) ,C2 (72.20) ,C3 (72,20) ,C4 (72,20) 
DO 20 J-l.N 
DO 20 1-1,72 

A 1 ( I , J) - ( 1 . 0 + 32 . 2095 *A ( I , J) ) / ( 1 . 0 + 52 . 85 *A ( I , J) ) 

A2 (I , J)- (0.534874+199. 0*A (I , J) -1990. 63*A (I , J) *A (I , J) )/ (1 .0+ 

1 333*244*A (I , J) ) 

A3 (I ,J)-(1.0+74.144*A(I ,J))/ (0.43368+24. 7442*A (I ,J)) 

20 CONTINUE 
DO 30 J- 1 .N 
DO 30 1-1,72 

WW(I,J)-(A1 ( I , J) +A2 ( I , J) *SQRT (C ( I , J) ) ) 

1 / (1 .0+A3 (I , J) *SQRT (C ( I , J) ) ) 

WW(I , J)-AMAX1 (WW(I ,J) ,0.0) 

WV(I ,J) -1.0/ (1.0+9. 2241 1*SQRT(B(I,J))+33.1236*B(I , J)+l 76.396 
1 *B (I , J) *B (I , J) ) 

WV(I ,J)-AMAX1 (WV(I , J) ,0.0) 

30 CONTINUE 
DO 40 J-l.N 
DO 40 1-1,72 

A (I , J)-AMIN1 (A ( I , J) ,0.06) 

B (I , J) -AMIN 1 (B (I ,J) ,2.0) 

C(l ,J)«AMIN1 (C(l ,J) ,8.0) 

40 CONTINUE 
DO 50 J-l.N 
DO 50 1-1,72 

YV ( I , J) -0 . 1 * (0 . 085 1 069*SQRT (B ( I , J) ) +0 . 32 3 1 05*B ( I . J) -0 . 1 87096 
1 *B ( I , J) *SQRT (B ( I , J) ) ) 

ZV ( I , J) -C . 00 1 * (0 . 239 1 86 *B ( I , J) - 0 . 092 2289*B ( I , J) *SQRT (B ( I , J) ) - 
1 0.0167413*8 (I , J) *B (I , J) ) 

50 CONTINUE 
DO 60 J— 1 , N 
DO 60 1-1,72 

B 1 (I ,J)- ( 5 . 6383 E- 4 + 1 .05173** (I ,J)- 39 . 0722 *A(I , J) *A (I , J) ) 

1 / (1 .0+202 .357*A (I , J) ) 

B2 ( I , J) - (0 .0779555+ 1 * • 40720*A ( I , J) + 3 . 1 585 1 *A ( I , J) *A (I , J) ) 

1 /(1.0+40.2298*A(I,J)) 


0032 B3 ( l , J) - (-0 . 038 1 305-3- 6368i»*A ( l , J) +7 • 9895 1 *A ( 1 , J) *A ( 1 , J) ) 

1 / (l . 0 + 62 . 5692 ** (l , J) ) 

0033 B4 (I , J) » (0. 0062 1 039+0 . 7 1006 1 *A ( I , J) -2 . 8524 1 *A ( I . J) *A ( I , J) ) 

1 / ( 1 . 0 + 70 . 2912 ** (I . J)) 

0031* YW ( I , J) -0 . 1*(B1 (I , J)+B2 (I , J) *SQRT (C ( I , J) )+B3 ( I . J) *C (I , J) 

1 +BMI,J)*C(I . J) *SQRT (C (I . J) ) ) 

0035 60 CONTINUE 

0036 DO 70 J-l.N 

0037 DO 70 1-1,72 

0038 C 1 ( I , J) « (-2 . 99542E-4+0 . 2382 1 9*A ( I , J) +0 . 5 19261»*A ( I . J) *A ( I , J) ) 

1 / ( 1 .0+10.7775*A (I , J) ) 

0039 C2 (I . J) - (-0.0291 325-2 . 30007*A ( I , J) +1 0 . 9460*A (I , J) *A (I , J) ) 

1 /(1.0+63.519*A(I . J) ) 

001*0 C3(l ,J)-(0. Oil* 3812+1. 80265*A (I ,J) -10.13H*A(I ,J)*A(I ,J)) 

1 /(1.0+98.1*758*A(I,J)) 

001*1 Cl* (I ,J)« (- 0 . 00239016 - 0 . 37 H 27 *A( I ,J)+ 2 . 35 M» 3 *A(I ,J)*A(I ,J)) 

1 / ( 1 . 0 + 120 . 228 *A (I , J) ) 

001*2 ZW(I ,J) -0.001* (Cl (I , J) +C2 ( I , J) *SQRT (C ( I , J) )+C3 (I , J) *C (I , J) 

1 +C1* (I , J) *C ( I , J) *$QRT (C ( I , J) ) ) 

001*3 70 CONTINUE 

001*1* RETURN 

001*5 END 


UNCOLLAPSIBLE LOOPS REASON FOR NON-VECTOR I ZAT I ON 


LINE 0018 
LINE 0017 
LINE 0011 
LINE 0010 


LINE 0021 
LINE 0018 
LINE 0015 
LINE 001 1 


STATEMENT CONTAINS NONVECTOR I ZABLE FUNCTION 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 
STATEMENT CONTAINS NONVECTOR I ZABLE FUNCTION 
NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 


NUMBER OF LOOPS IN ROUTINE - 0012 
NUMBER OF VECTOR I ZABLE LOOPS = 0008 


NUMBER OF STACKLIBED LOOPS - 0000 


0001 SUBROUTINE TRC03 (N,A,B,C,D,E,F,TAUF ,TXUF) 

C 

C 

0002 DIMENSION A (72 , 20) , B (72 , 20) ,C (72,20) ,D (72,20) ,E (72,20) ,F (72,20) 

0003 DIMENSION TAUF (72,20) ,TXUF (72,20) 

0001. DIMENSION XV(72,20) ,XW(72,20) ,XZ (72,20) ,DX (72,20) , DY (72,20) 

0005 DO 10 J-l.N 

0006 DO 10 1=1,72 

0007 XV ( I , J) -3 . 1 *A ( I , J) / ( 1 .+ 1 5 . 1 *A ( I , J) **0 . 56) 

0008 XW(I ,J)-0.01**B(I ,J) / (1 .+0 . 9*B ( I , J) **0.57) 

0009 DX ( I , J) - 6 . 7*C ( I , J) / ( 1 .+16 . *C ( I , J) **0 . 6 ) 

0010 DY (I , J) -27»*D (I , J) **0.83 

0011 XV (I , J) -XV ( I , j) +DX ( I , J) +DY ( I , J) 

0012 XW ( I , J) -XW ( I , J) +DX ( I , J) +DY (I ,J) 

0013 XV ( I , J) -AMIN 1 (XV ( I ,J) ,30.) 

0011* XW(I ,J)-AMIN1 (XW ( I ,J) ,30.) 

0015 TAUF (I , J) -O. 38 L 61 5381**EXP (-XV (I , J) ) + 0 . 6 15381*6 15 *EXP (-XW(I ,J)) 

0016 DY(I,J)-9.79*D(I,J) 

0017 DY (I ,J)-AMIN1(DY (I , J) , 30 .) 

0018 F (I ,J)-F (I ,J) / (E (I ,J) *1013.25) 



0019 

XZ ( 1 , J) -1* . *3A5 . 28*E ( 1 , J) / (0 . 879&*F ( 1 , J) ) 


0020 

XZ ( 1 , J) —1* . 398*F ( 1 , J) * (SQRT (1 .+XZ (1 , J) ) -1 .) 


0021 

XZ (i , J) > 1 . 0 - 81 . 2 1 / 120 . 0 *( 1 . 0 -EXP (XZ ( 1 ,J))) 


0022 

TXUF ( 1 . J) -XZ ( 1 » J) *EXP (-DY ( 1 , J) ) 


0023 

10 CONTINUE 


0021 * 

RETURN 


0025 

END 


UNCOLLAPSIBLE LOOPS REASON FOR NON-VECTORIZATION 


LINE 0006 

LINE 0017 STATEMENT CONTAINS NONVECTOR 1 ZABLE 

FUNCTION 

LINE 0005 

LINE 0006 NONVECTOR 1 ZABLE LOOP NESTED WITHIN 

LOOP 

NUMBER OF 

LOOPS IN ROUTINE = 0002 


NUMBER OF 

VECTOR 1 ZABLE LOOPS * 0000 


NUMBER OF 

STACKLIBED LOOPS * 0000 



0001 

C 

SUBROUTINE CLOUDS 

0002 

C 

COMMON/CSC/CSS (72,20) ,CCU(72,20) ,CH (72,20) ,CC (72,20) 
1, CS (72,20,20) 

0003 


COMMON/ 1 N 1 T/NP , NPP 1 , NPM 1 

0001 * 


DIMENSION CT (72,20) ,CU(72, 20) ,CR (72,20,20) 

0005 


DO 10 IP-1, NP 

0006 


DO 10 IL-1,72 

0007 


CH (IL, IP)=1 .0 

0008 


CC (1 L, IP)-1 .0 

0009 


CT(IL. IP)-1 .0 

0010 


CU ( 1 L , 1 P) *1 .0 

0011 

10 

CONTINUE 

0012 


DO 15 IP-1, NP 

0013 


DO 15 IX-l.NP 

OOlL 


DO 15 IL-1,72 

0015 


CS ( 1 L, 1 X, 1 P) -1 .0 

0016 


CR (IL, IX, 1 P) -1 .0 

0017 

15 

CONTINUE 

0018 


DO 20 1 P- 2 , NP 

0019 


DO 20 IL-1,72 

0020 


CT ( 1 L, 1 P) -1 . -CCU ( 1 L, IP-1) 

0021 


CT ( 1 L , IP) -AMIN 1 (CT ( 1 L , 1 P- 1) , CT ( 1 L , 1 P) ) 

0022 


CH ( 1 L , 1 P) -CH ( 1 L , 1 P- 1 ) * ( 1 .-CSS (IL, IP-1)) 

0023 

20 

CONTINUE 

0021 * 


DO 30 IP- 2 .NP 

0025 


DO 30 IL-1,72 

0026 


CH (1 L, IP) -CH (1 L, IP) *CT(IL, IP) 

0027 

30 

CONTINUE 

0028 


DO 1*0 l-l.NPMl 

0029 


IP-NP-I 

0030 


DO 1*0 IL-1,72 

0031 


CU ( 1 L , IP )«1 ,-CCUdL, IP) 

0032 


CU (1 L, 1 P) -AMI N1 (CU (1 L , IP+1) ,CU (1 L, 1 P) ) 

0033 


CC(IL.IP) -CC (J L , 1 P+l) * (1 .-CSS(IL, IP)) 

0031 * 

1*0 

CONTINUE 



0035 


DO 50 l-l.NPMl 

0036 


1 P-NP- 1 

0037 


DO 50 IL-1,72 

0038 


CC(IL.IP) -CC(IL.IP) *CU (1 L, IP) 

0039 

50 

CONTINUE 

0060 


DO 60 IP-2.NP 

0061 


IPM1-IP-1 

0062 


DO 60 IX-1.IPM1 

0063 


DO 60 IL-1,72 

0066 


CR ( 1 L, 1 X, 1 P) *1 .-CCU (1 L, 1 P- 1) 

0065 


CR (It, IX, IP) -AMIN1 (CR (1 L, 1 X, 1 P-1) ,CR (IL, IX, IP)) 

0066 


CS (IL. IX. IP) -CS(IL, IX,IP-1)*(1.-CSS(II.,IP-1)) 

0067 

60 

CONTINUE 

0068 


DO 70 IP-2.NP 

0069 


IPM1-IP-1 

0050 


DO 70 IX-1.IPM1 

0051 


DO 70 IL-1,72 

0052 


CS (1 L, IX, 1 P) -CS (IL, IX, IP) *CR (IL, IX, IP) 

0053 

70 

CONTINUE 

0056 


DO 80 I-1.NPM1 

0055 


1 P-NP- 1 

0056 


IPP1-IP+1 

0057 


DO 80 IX-IPPl.NP 

0058 


DO 80 IL-1,72 

0059 


CR (1 L, 1 X, 1 P) -1 .-CCU (1 L, IP) 

0060 


CR(IL, IX, IP) —Art 1 N 1 (CR ( 1 L, 1 X, 1 P+1) ,CR (I L , 1 X, IP) ) 

0061 


CS (1 L, 1 X, 1 P) -CS (1 L, 1 X, 1 P+1) * (1 .-CSS (1 L, IP)) 

0062 

80 

CONTINUE 

0063 


DO 90 l-l.NPMl 

0066 


1 P-NP- 1 

0065 


IPP1-IP+1 

0066 


DO 90 IX-IPPl.NP 

0067 


DO 90 IL-1,72 

0068 


CS(IL, IX, IP)«CS(IL, IX, IP)*CR(IL, IX, IP) 

0069 

90 

CONTINUE 

0070 


RETURN 

0071 


END 


UNCOLLAPSIBLE LOOPS 

REASON FOR 

NON-VECTOR IZAT ION 

LINE 0063 

LINE 

0066 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0056 

LINE 

0057 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0057 

LINE 

0061 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0058 

LINE 

0060 

STATEMENT CONTAINS NONVECTOR 1 ZABLE FUNCTION 

LINE 0068 

LINE 

0050 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0060 

LINE 

O 

O 

ro 

LOOP WITH VARIABLE INITIAL/TERMINAL VALUE 
NESTED WITHIN LOOP 

LINE 0063 

LINE 

0066 

RHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 

LINE 0062 

LINE 

0063 

NONVECTOR 1 ZABLE LOOP NESTED WITHIN LOOP 

LINE 0035 

LINE 

0038 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0028 

LINE 

0033 

PROPERTY OF EMBEDDED LOOP PROHIBITS 
VECTOR IZAT ION OF LOOP 

LINE 0030 

LINE 

0032 

STATEMENT CONTAINS NONVECTOR 1 ZABLE FUNCTION 

LINE 0019 

LINE 

0022 

RHS ARRAY HAS POSSIBLY RECURSIVE PROPERTIES 



LINE 0018 
LINE 0012 


LINE 0019 NONVECTOR I ZABLE LOOP NESTED WITHIN LOOP 
LINE 0013 LOOP WITH VARIABLE INITI AL/TERHINAL VALUE 
NESTED WITHIN LOOP 

NUMBER OF LOOPS IN ROUTINE - 0025 
NUMBER OF VECTORIZABLE LOOPS - 0011 
NUMBER OF STACKLIBED LOOPS - 0000 


51 



